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I. INTRODUCTION 



A. BACKGROUND 

For more than two decades, the need to solve control 
problems in aerospace application has been the primary 
driving force behind the modern control theory development. 
Problems in manoeuvring, guidance and tracking of aircrafts 
and space vehicles have motivated the development of various 
control design and synthesis methods. One of the methods, 
the so-called Linear Quadratic Control or LQ-control 
[Refs. 1,2] has been widely used and is treated extensively 
in the control literatures. Unlike most classical methods 
where the design are based on conventional time response 
criteria, the LQ theory treats the problem of designing 
controllers as that of minimizing a quadratic cost function 
of states and control inputs. The design problems become 
that of selecting suitable weighting matrices in the 
performance index. Two questions naturally arise from this 
method. First, how does one select the weighting matrices 
and second, once a set weighting matrices is selected, how 
does one know that it is a good design. There are generally 
two approaches to the first problem; the obvious one is to 
rely on physical arguments and a certain amount of trial and 
error [Refs. 3,4] Unfortunately, such formulation can be 
obtained in only a few cases . Reference 4 provides a few 
guidelines that can be employed. The second approach is to 
avoid the physical aspect of the performance index but 
instead try to relate the weighting matrices with some other 
performance specification. For example, Tyler and Tuteur 
[Ref. 5] expressed the characteristic polynomial as an 
explicit function of the weighting matrices for single input 
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single output (SISO) system with diagonal weighting 
elements. Root-Locus type of procedure were used to show 
how the variation in weighting elements affect the eigenva- 
lues of the closed-loop system. Similar relations were 
explored in [Refs. 6,7] in which more general expressions 
were obtained. Their uses were restricted to single input 
case due to difficulty in handling polynomial matrices. 
Solheim [Ref. 8] later developed a sequential design proce- 
dure based on diagonalized (decoupled) system. More recent 
results of eigenvalue placement in optimal control problem 
are presented in [Refs. 9,10,11,12]. 

The second issue of LQ design, i.e., whether a good 
design has been obtained once the performance index has been 
fixed, is related to the multivariable nature of system. 
Unlike the single input case where the closed- loop eigenva- 
lues uniquely define the feedback gain and hence the 
weighting matrices, the MIMO structure provides additional 
dimension which allows further tradeoff for properties other 
than the closed-loop eigenvalue location. An example is the 
gain and phase margin, or in MIMO case, the so-called 
robustness criteria. Robust-control has been the subject of 
extensive researches in recent years and results relating to 
optimal control can be found in [Refs. 13,14,15]. 

It becomes evident from the above discussion that LQ 
control design and synthesis are not just a matter of speci- 
fying performance index. Optimal in the sense of satisfying 
performance index and perhaps pole location does not neces- 
sary means that a good design has been obtained. Other 
criteria like disturbance rejection, robustness and sensi- 
tivity need to be considered and incorporated in the design 
procedure. This is the motivation behind our present 
research which in turn leads to the development of the 
synthesis package. An overview of the thesis is described 
in the next section. 
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B. OVERVIEW 



In this section an overview of the thesis is given. A 
background of multivariable optimal control theory in terms 
of its structure, frequency domain characteristic and asymp- 
totic properties is first discussed together with derivation 
of some useful relationships in Chapter 2. General robust- 
ness concepts and its application in Linear Quadratic (LQ) 
Control is presented in Chapter 3. In Chapter 4, a computer 
aided design package for pole-placement synthesis based on 
numerical optimization technique is presented. This package 
provides a useful computational tool to support the material 
in the remaining chapter. A step-by-step pole-placement 
synthesis procedure is also presented to illustrate how the 
package can be used to design optimal control system that 
meet time response criteria and other properties. The use of 
the pole-placement synthesis package in two actual design 
problems is demonstrated in Chapter 5. Results in terms of 
frequency and time response properties , robustness (singular 
value decomposition) are compared to those obtained by other 
design methods. Program listings and an example of a design 
session are included in the Appendix. 
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II. MULTIVARIABLE LINEAR QUADRAT IC CONTROL 



In this Chapter a brief review of Multivariable Linear 
Quadratic Control theory is presented. Structure of LQ 
system is first given, general stability properties are then 
discussed. In section two the general pole assignment 
problem is formulated for the MIMO state feedback system. 
Interpretation of the state and control input weighting 
matrices and their effect on the closed-loop time response 
behaviors and the asymptotic properties are presented. 

A. LINEAR QUADRATIC SYSTEM 

Consider the following linear time-invariant state-space 
system given by the equations: 

x(t) = Ax(t) + Bu(t) (eqn 2.1) 

y ( t ) = Cx ( t ) (eqn 2.2) 

where x(t) is the n-diinensional state vector, u(t) is the 
m- dimens ional control vector and y(t) is the p- dimens ional 
output vector. A ,B and C are real constant matrices of 
dimension nxn , nxm and pxn respectively. Assuming that they * 
form a controllable pair(A,B) and an observable pair(A,C), 
the optimal feedback control law is obtained by minimizing 
the following quadratic performance index. 

f 00 

J = ) ( x T Q x + p u T R u)dt (eqn 2.3) 

^0 

where R is positive definite (R>0) for bounded input and p 
is a scalar. Q is a semi-positive definite matrix (Q^O). 
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When both Q and R are diagonal matrices, P defines the 
relative weight between the state and control weighting 
matrices in the performance index. The quadratic form x T Qx 
and u T Ru provide a weighted measure of the magnitude of the 
states and control vector respectively. The steady state 
control law that minimizes J is given by. 



u ( t ) = -Fx(t) 



(eqn 2.4) 



where F is the feedback gain matrix which is given by 

F = -R“Vp (eqn 2.5) 

The positive definite matrix P is given by the solution to 
the steady state Riccati equation. 

PA + A t P + Q - PBR'Vp = 0 (eqn 2.6) 

Equation 2.5 and 2.6 are well-known results in optimal 
control theory that yield the optimal closed-loop system 

x ( t ) = [A + BF ] x ( t ) (eqn 2.7) 

whose closed-loop poles are given by the eigenvalues of the 
matrix [A+BF] . The LQ system given by equations 2.1 through 
2.6 is closed-loop stable and can be represented by the 
feedback configuration in both the time and frequency domain 
as shown in Figure 2.1 . 
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Frequency Domain 




Time Domain 



Figure 2.1 State Feedback System-Time and Frequency Domain. 
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In Figure 2.1, if the loop is broken at the input (X) as 
shown, the loop transfer function is given by 



G (s ) = F ( si 



A) _1 B 



( eqn 2.8) 



The matrix [I + G(s)] is called the Return Difference Matrix 
and will be shown to have some important feedback properties 
in the following section. 



B. 



POLE ASSIGNMENT PROBLEM 



In the most general form, the state feedback pole 
assignment problem in control system design can be formu- 
lated precisely as follows: 

"Given real matrices (A,B) of order (nxn, nxm) respec- 
tively and a set of n complex numbers ( A, , A, , . . . A n ) closed 
under complex conjugation. Find a real mxn matrix F such 
that the eigenvalues of [A + BF] are A ; , i = 1, 2,3... n." 

In general, the closed-loop eigenvalues of (A+BF) can be 
arbitrary located in the complex plane, with the only 
restriction that complex characteris tic eigenvalues must 
occur in complex conjugate pair. In other words, if the 
matrices (A,B) are a completely controllable pair, the 
stability of the system can always be improved by state 
feedback. If the (A,B) pair is not completely controllable, 
then the system is required to be ' stabilizeable ' meaning 
that in its controlability canonical form given by equation 
2.9 below, the matrix is asymptotically stable or any 
unstable subspace of equation 2.9 is also in the 
controllable subspace. 
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X, (t ) 




A H 


A,, 


X x (t ) 




1 

o 


An 



x , (t ) 


+ 


B 




u,(t) 


x 2 (t)_ 




0 




u*(t) 



( eqn 2.9) 

The solution to the pole assignment problem in the 
single input (m=l) case, when it exits, can be shown to be 
unique. In the multiple input case (l<m<n), the solution of 
the so-called inverse eigenvalue problem is, in general, 
undetermined with many degrees of freedom. Additional 
conditions must be supplied in order to eliminate the extra 
degrees of freedom. This has been an area of active 
research in recent years and a number of approaches have 
been developed to relate the extra degrees of freedom with 
properties such as system eigenvectors, transmission zeros 
and robustness [Refs. 16,17]. In this work, it is shown 
that Linear Quadratic formulation incorporating equation 2.3 
will partly remove the uncertainty that exists in the 
multiple inputs case. It will be shown that useful proper- 
ties like robustness are guaranteed . It is also shown that 
the LQ type of pole placement formulation, when combined 
with eigenvectors assignment type of formulation, will 
produce some very useful design and synthesis procedures. 
In the present work, however, only the LQ eigenvalue place- 
ment problem is addressed; the problem of combined eigen- 
value and eigenvector assignment is briefly described in 
Chapter 5. 
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C. WEIGHTING MATRICES AND SYSTEM PERFORMANCES 



This section briefly reviews the effect of weighting 
matrices on system performance. The physical aspects of the 
weighting matrix for both single and multiple input systems 
are presented first together with some discussion on the 
asymptotic behavior of the LQ system. 

It was shown in the last section that under complete 
controllable conditions, the t ime- invariant linear system 
can be stabilized by a linear feedback control law. For the 
regulator- type problem where the aim is to bring the system 
from an arbitrary initial state to the zero state, the 
closed loop poles can be chosen far to the left on the 
complex plane. Convergence to the zero state is fast but 
the large input required may not be practical. This natu- 
rally leads to an optimization problem where the trade off 
is between speed of convergence to zero and the magnitude of 
the input. This is reflected in the two quadratic terms in 
the performance index. The quantity x r Qx in the first term 
of the performance index is a measure of the extent in which 
the state at time t deviates from the zero state. The 
matrix Q determines how much weight is attached to each of 
the component of the state. The integral u /'(x T Qx)dt is a 
criteria for the cumulative deviation of x(t) from the zero 
state during the interval. 

The problem of large control input is resolved by 
incooperating the second quadratic term, J"(u T Ru)dt. Larger 
value in the element of the control weighting matrix R will 
result in smaller input. It will be shown later that one can 
manipulate R to achieve some secondary design objectives. 
The remaining parameter p , that need to be specified, 
accounts for the relative weighting between the state and 
the control inputs. Selecting optimal value for p depends 
on the particular problem and the design requirement. As p 
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decreases, the integrated square regulating error decreases 
but the integrated square input increases. Very often, the 
optimal control problem is solved for many different values 
of P A graph similar to that shown in figure 2.2 is 
obtained where the integrated square regulating error is 
plotted versus the integrated square input. An appropriate 
value of p is then selected to give sufficiently small 
regulator error without excessively large input. 




Figure 2.2 Selecting Relative Weighting Parameter p . 
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The special case where p approaches zero ,the so-called 
asymptotic properties, has been shown to provide good 
insight for LQ control system design. Some of the results 
from [Refs. 4,18] are summarized below. 

As p decreases to zero for the system given by equa- 
tions 2.1 through 2.6, some ( say q ) of the closed- loop 
poles go to infinity while other ( n-q ) stay finite. Those 
remaining finite approach the left half plane zero of 
det [B T ( - s I - A T ) Q (sI-A) 1 B] . If m is the dimension of the 
input vector, then at least m closed-loop poles approaches 
infinity. All closed-loop poles that go to infinity do so 
by grouping into several Butterworth patterns of different 
order and radii. For p approaching infinity, the closed- 
loop poles approach the mirror image of the plant open- loop 
poles. To illustrate the asymptotic concept, two examples 
from [Ref. 4] are given below. 

An example of the asymptotic properties of a single 
input system is shown in Figure 2.3 . The closed-loop poles 
of a position control system using LQ type feedback are 
plotted as a function of p . The system has two open- loop 
poles at -4.6 and 0.0 . As p decreases, the closed-loop 
poles go to infinity along two straight lines that make an 
angle of 7r /4 with the negative real axis. As p 

approaches infinity, both closed-loop poles first meet on 
the negative real axis and then approach the open- loop 
poles at (-4.6, 0 ) . 

Figure 2.4 shows the asymptotic loci of the closed-loop 
poles for a multiple input system where n=4 and m=2. There 
are four open-loop poles at (-0.006123, ±j0. 09353) and 
(-1.250, ±jl.394) and both the state and control weighting 

matrices are taken to be of a diagonal form. As p 
approaches to zero, one closed-loop pole stay finite at 
-1.002. The remaining three go to infinity, two of which 
assume a second order Butterworth pattern and the last one 
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approaches on the negative real axis. When p approaches 
infinity, all closed- loop poles approach the open-loop 
poles . 

Many researchers have explored this asymptotic proper- 
ties. Stein [Refs. 9,10] developed a procedure to select 
the control weighting matrices for a desired closed-loop 
asymptotic eigenstructure . It will be shown in the subseq- 
uent Chapters that the asymptotic properties provide useful 
guideline for the design procedures to be described. 




Figure 2.3 Single Input Asymptotic Root Loci. 
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Nearby Poles 



Figure 2.4 Multiple Input Asymptotic Root Loci. 
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III. ROBUSTNESS THEORY AND LINEAR QUADRATIC CONTROL 

In the last Chapter, the effect of weighting matrices on 
the closed-loop pole of the LQ system was discussed. This 
Chapter addresses yet another important feedback property 
that control system designers are concerned with: 
Robustness. Robustness theory was developed when it was 
realized that the classical single loop Nyquist test was not 
adequate to guarantee stability when the multivariable open 
loop plant deviates from its model due to a variety of 
reasons [Ref. 15]. In the following section , the concept 
of robust design for both SISO and MIMO general feedback 
system is reviewed. The singular value analysis is 
discussed in term of multiplicative type of disturbances. 
Finally, robustness for linear quadratic state feedback are 
presented in terms of the effect of weighting matrices on 
the singular value curve. 

A. ROBUSTNESS CONCEPTS 

The concept of robustness for SISO system can best be 
described in terms of the definition of phase and gain 
margin. A system characterized by good gain and phase 
margin implies that changes in the plant model parameters 
and changes in the loop gain and/or phase may be accommo- 
dated without loss of stability. The gain and phase margins 
of a SISO system are defined with reference to the perturbed 
system in Figure 3.1 

Assuming that the unperturbed system (l(jw)=l) is 
stable, the positive (or negative) phase margin is the value 
of <t> greater (or less) than zero at which the perturbed 
system with l(jw) = k*exp(j$) becomes unstable. The upward 
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I 




(or downward) gain margin is the smallest (or greatest) 
value of 1 ( j w) =k= cons tant , for which the system become 
unstable. With reference to the classical Nyquist plot in 
Figure 3.2 , gain and phase margin are defined as, 

GM f (upward) = l/k ( GM l (downward) = l/k 2 

PM + = a , PM ' = a z 
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A set of minimum guaranteed GM and PM may be obtained by 
defining “ 0 = min[l + g(jw)] as shown in Figure 3.3 where, 

GM = 1/(1 ± « 0 ) 



and 



- 1 , 



PM = ±cos ( 1- 



a o/2) 



Imcuj 




Figure 3.2 The Nyquist Plot - Gain and Phase Margin. 
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Note that in the above definition of GM and PM, either 
<!>! 'or' Kj is allowed to change in the loop. The allowed 
changes are therefore very restrictive. A more useful defi- 
nition of the gain and phase margin for MIMO system that 
accounts for simultaneous changes in both ^ and/or K ; (the 
so-called universal gain and phase margin) has been derived 
in [Ref. 19], From Figure 3.3 it can be seen that <* 0 may 
provide a rather conservative estimate of the actual gain 
and phase margin. 




Figure 3.3 The Minimum Guaranteed GM and PM. 
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In MIMO system, gain and phase margin characterize the 
ability of the system to tolerate gain and/or phase changes 
within all loops simultaneously. Figure 3.4 shows a 

perturbed MIMO system with the assumption that L(jw) = Dia(l, 

( jw) , 1-2. ( jw) . . . ) . As in the SISO case, the system will 
remain stable as longer as 1; (jw) satisfy GM | < lc ; < GMf 

(assuming that 4> ; s are zero ) . For the case when the 

magnitude of l(jw) are constant, the system will remain 

stable for all 4> ; s that satisfy PM”< < PM*. 




Figure 3.4 The Perturbed MIMO System. 
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The set of minimum guaranteed multivariable gain and 
phase margins are determined by the minimum singular value 
( £_ ) of the return difference matrix [I + G(s)], where G(s) 
is the open loop transfer matrix and consists of the plant 
and its controller. Note the similarity between the SISO 
and MIMO cases, 



a o : nearness of (1 + g(s)) to the origin. 

: nearness of matrix [I + G(s)] to singularity. 

Two important results that are related to multivariable 
phase and gain margin developed in [Ref. 15] are next 
presented as a theorem. 

THEOREM: The multiplicative perturbed system (Figure 

3.4) is stable if either of the following conditions hold: 



1 . 


1 [I ♦ G(s)] > 


- [L" (s) 


- I] 


(eqn 


3.1) 


2. 


£[t*G''(s)] > 


a [L (s) - 


I] 


(eqn 


3.2) 



where a , a_ denote the maximum and minimum singular values 
of [I+G(s)] respectively. 



It will be shown that condition 1 can be related to the 
optimality condition of LQ system and hence provides a 
useful relationship between the weighting matrices and 
robustness . 
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B. ROBUSTNESS IN LQ SYSTEM 



Robustness properties pertaining to LQ system are 
closely related to the frequency domain optimality condi- 
tion. For the SISO case, Kalman [Ref. 1] showed that the 
return difference transfer function satisfies the inequality 

( 1 + g(jw) ) > 1 ( « 0 = 1) (eqn 3.3) 

Inspection of the Nyquist diagram in Figure 3.5 clearly 
indicates that a SISO LQ state feedback has a guaranteed 
infinite upward gain margin , 0.5 downward gain margin and a 
minimum phase margin of ± 60 deg. 




Figure 3.5 SISO LQ Nyquist. 
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In a similar manner, multivariable stability margin can 
be derived from the multivariable version of condition 1 in 
equation 3.1 ; 

[I + G(s)] H R [I + G ( s ) ] > R (eqn 3.4) 

where 

G (s ) = R' 1 B T P[sI - A] 1 B 

and P > 0 satisfies the steady state Riccati equation. 

It can be shown that condition 1 in equation 3.1 above 
can be written as 

Z. [I + R^ 2 G(jw) R ,/Z ] > 1 (eqn 3.5) 

which can be reduced to 

£_[I + G ( jw) ] > 1 (eqn 3.6) 

where R is diagonal ( i.e. R = p I for some positive 
scalar p). Like the SISO case, the multivariable LQ regu- 

lator with loop transfer matrix G(s) that satisfies equation 
3.5 and 3.6 has, in all the feedback loops, a guaranteed 
minimum gain and phase margin given by, 

GM = 1/2, and PM = +60 deg 

The case when R is not diagonal is interesting as 
condition 3.6 no longer apply and additional trade-off can 
be obtained by including off diagonal terms in the control 
input weighting matrix. This is especially useful when the 
nature and structure of the disturbances are known. Design 
involving the off diagonal R matrix will be shown in 
Appendix A. 



31 



IV. DESIGN PROCEDURE 



A. GENERAL 

The background material described in Chapter 2 and 3 are 
the basis for developing a computer aided design package and 
the corresponding design procedure. The design philosophy 
presented in this Chapter is most useful when a designer is 
able to characterize the desired system in terms of closed- 
loop eigenvalues and time response. Some initial physical 
insight of the weighting matrices , their asymptotic proper- 
ties and the nature of the perturbation will be useful in 
the design process. This Chapter begins with a discussion 
of the various approaches to the LQ pole placement problem. 
The selected approach and the corresponding computer aided 
design package is presented. The design philosophy and 
design procedure for both the reduced order and full order 
model are then given. 

B. APPROACHES 

All the pole placement algorithms for LQ system that 
have been developed so far require an expression that relate 
the characteristic equation of the optimal system with the 
elements of the weighting matrices. Two such formulations 
for the stabilizeable and detectable t ime - invariant linear 
system (equation 2.1) and the quadratic criteria (equation 
2.3) are given by 

^c( s )^( -s ) = $ ( s ) <t> ( - s )det [I + R - ’ H T (-s) Q H(s)] (eqn 4.1) 
and 
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.X T - 1 

<£ c (s)<^(-s) = 4>(s )<f>( -s )det [1+ 1/pR 2 H ( -s )QH(s )R 2 ] (eqn 4.2) 

where <£ c (s) = det[sl -A +BF] and </>( s) = det[sl - A] are 

the closed-loop and open- loop characteristic polynomials 
respectively. H(s) = C(sl-A) - ^ B is the open loop transfer 
matrix of the system. 

Both formulations have been used in root-locus type of 
procedure to investigate how the closed- loop poles move as 
weighting matrices changes [Refs. 5,6]. For MIMO case, 
there has been little success due to problems involving 
polynomial matrices. In this work, a different approach is 
adopted. Equation 4.1 or 4.2 is formulated as a numerical 
optimization problem in which the objective function is made 
equal to the determinant part of equation 4.1 or 4.2 

Obj = det [I ♦ R-'h T (s) Q H(s)] (eqn 4.3) 



or 



_x T _JL 

Obj = det [I + l/ p R~ 2 H T (-s) Q H(s)R Z ] 



( eqn 4.4) 



For a given desired closed loop pole s = s^ , equation 
4.3 or 4.4 becomes 



^ c ( s d ) ~ s d ) = 0 (eqn 4.5) 

Providing that <t>(s d )</>(- s d ) is not equal to zero, the 
objective function must equal to zero if the particular Q, R 

I 

set is to correspond to the desired closed-loop poles. 
Convergence to zero for a given set of Q and R is therefor 
automatically guaranteed. 

The pole placement problem can therefore be solved as an 
unconstrained multivariable optimization in which the 
elements of Q and R are varied to make the objective 
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function in equation 4.3 and 4.4 approach zero. This was 
done during the early phase of the work. It was later 
discovered that more insight to the problem can be obtained 
by first transforming the problem to an appropriate coordi- 
nate system and then to perform pole placement one at a time 
This is of advantage as a system designer is often satis- 
fied with several open-loop poles in a large system. 
Reassigning poles in the the reduced- order model will reduce 
computer time and memory requirement. 

The optimization routine selected for this work is the 
so-called SUMT method (Sequential Unconstrainted 
Minimization Techniques) obtained from the ADS package in 
[Refs. 20,21]. In this method, the objective function 

(eg. equation 4.3 or 4.4 ) and any constraint equation are 
formulated into an augmented objective function in which the 
problem is solved as an unconstraint optimization task. 

C. POLE PLACEMENT ALGORITHMS 

For ease of implementation and better insight, only the 
problem of determining the state weighting matrix Q (given 
R) that gives a set of closed-loop eigenvalues is consid- 
ered. It must be emphasised that the algorithm can also be 
formulated to determine R (for a given Q) or to vary Q and R 
at the same time. In most cases, the present formulation is 
adequate as designers usually have some knowledge about the 
control weighting matrix. System matrix that has real and 
distinct eigenvalues is presented first, follows by cases 
where A has complex eigenvalues and repeated eigenvalues. 

1. S ystem with Real and Dis tinct Eig envalues 

The original system given by equation 2.1 and 2.2 is 
first transformed into a diagonal form using the 
transformation given by; 
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x(t) = Mz(t ) 



(eqn 4.6) 



where M is an eigenvector matrix corresponding to the system 
matrix A. and z(t) is the new state vector. The trans- 
formed system in the new coordinate is given by, 

z(t)= Az(t)+MBu(t) (eqn 4.7) 

where A is a diagonal matrix dia[ A, , A, , X 3 , . . . ] 

The performance index (equation 2.3), when expressed 
in terms of the new state vector z(t) becomes, 



J = J Q ( x t Qx + u T R u ) dt 
= J q (z T M t QM z + u T Ru) dt 

= \ ( z r Qz + u T R u ) dt (eqn 4.8) 

0 

A T 

where Q = M QM. 

It can be shown that to move an open- loop pole to 

/\ 

its new location given by s ; = s d , only Q; is required and 
other Qs ' have no effect on the pole assignment. As an 

example, to move the open- loop pole at s = X ^ for 
A = dia[ X, , X x , Xj , . . .] to its new location s= \ n , only 

Q= dia [0 , q , 0, . . . ] is required. 

Q can then be selected according to equation 4.3 or 
4.4 using the optimization routine. As currently imple- 
mented in the program, there is no constraint equation 

A 

formulation. Once the value of Q that satisfies the desired 
pole location is obtained, the system is transformed back to 
the original coordinate system via, 



-Ta -I 

Q = M QM 



(eqn 4.9) 
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With Q known, and R given, the optimal feedback gain F can 
be obtained by solving the steady state Riccati equation as 
given in equation 2.5 and 2.6 . Since the pole placement is 
done in the decoupled coordinate system, only the eigenva- 
lues that correspond to Q ; is reassigned; all other eigenva- 
lues remain unchange. It can be also shown that the 
eigenvector corresponding to an eigenvalue is also unchange, 
this property will be shown to be useful in the reduced 
order formulation of the linear quadratic problem. 

If desired, the problem here can also be formulated 
to move more than one eigenvalues in one run. This can be 
done by modifying the objective function to include more 
terms as follows; 

Obj = 2 det [I + R' H ' (-s )QH(s )R* ] (eqn 4.10) 

i -i 

where n is the number of poles to be reassigned. 

The augmented matrix A 3ug = [A + BF] is then 
computed. If desired, the procedure may be repeated to move 
other open loop eigenvalue to its specified position using 
the new A aog as the starting plant matrix. This will in 
turn result in another set of Q and F . The effective Q e 

and F e after n reassignments are given by 

Q e = Q; + Q*. + • • • Q n ( eC l n *.11) 



and 



F e = F, + F t + . . . F n (eqn 4.12) 

The above pole assignment procedure can also be applied to 
an optimal system where an initial starting Q and R are 
given. A good example is when the control system designer 
has some knowledge of the weighting matrix but would also 
like to meet a specific time response requirement. 
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2 . System with Complex Eigenvalues 

If the same similarity transformation mentioned in 
the last section is used, the transformation matrix will be 
complex. To be able to work with real matrix, an auxiliary 
transformation of the form given by equation 4.13 is used; 



T = ML and M is the eigenvector matrix (equation 4.6) The 
matrix L is given by. 



x(t) = Tz(t) 



(eqn 4.13) 



0.5 - 0 . 5 j 0.0. 

0.5 0 . 5 j 0.0 



0.0 



L = 0.0 0.0 1.0 



1.0 



0.0 



1.0 



( eqn 4.14) 



The transformed system is then given by, 



z(t) = A Z(t) + TBu(t) 



( eqn 4.15) 



with performance index , 




(eqn 4.16) 



T 

where Q is now given by Q = T Q T. 
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It can be shown that to move a pair of complex 
eigenvalues given by s = a + b j , a weighting matrix of the 
form 



A 

Q = 



0 

0 

0 

0 

0 



0 

A 

0 



0 . 

0 

A 

V. 



0 



( eqn 4.17) 



with q> = q iH is required. 

A 

In a similar manner, Q can be obtained by using the 
optimization routine, with the condition q- = q formulated 
as a constrained equation. Inverse trails format ion and 

determination of Q e and F e are identical to the distinct 
eigenvalue case with M in equation 4.9 replaced by T. 



3. S ystem With Repeat ed Eigenval ues 

In this case, the system matrix cannot be diagonal- 
ized but the general procedure given above still apply. The 
system is first transformed into the Jordan canonical form ; 



J = U’*A U 



( eqn 4.18) 



where U is a transformation matrix which is not the eigen- 
vector matrix M. An example of the Jordan form is given 
below, 



A 



J = 



0 

0 



1 

A 

0 



0 

1 

A 



( eqn 4.19) 
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The only difference with the two procedures mentioned above 
is that the pole reassignment has to begin at the bottom of 
each Jordan block. For example, in the system given above 
(equation 4.19), the first re-assignment will result in a 
new system given by equation 4.20 



X 



1 



J = 



0 

0 



A 



0 




(eqn 4.20) 

Determination of Q and F using the optimization and the 
inverse transformation routines are identical to the 
distinct eigenvalues case with M in equation 4.9 replaced by 
U. 



D. COMPUTER PROGRAM DESCRIPTION 

The computer aided design package developed in this 
thesis is illustrated in Figure 4.1 The program is devel- 
oped using the top-down approach with special purpose 
subroutines called by the main control program. The driver 
program supports three independent modes with the data entry 
portion common to all modes: 

1. Data Entry : The system matrices A, B, C and/or F, Q, 
R etc are entered through a data file. Design vari- 
ables such as desired poles locations, elements of 
matrices to be varied etc, are also specified through 
the input data file. 

2. Pole Placement Mode : In this mode, an arbitrary set 

of closed-loop eigenvalues is assigned by selecting 
the appropriate state weighting matrix. As shown in 
Figure 4.1, the transformation matrices for various 
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cases are computed first. The pole placement is 
achieved using the numerical optimization routine 
described in the last section; Q is obtained and then 
inverse- trans formed to the original co-ordinate 
system. If desired, the results can be used as an 
input to the Linear Quadratic Control Program. 

3. Linear Quadratic Control Mode : This part of the 

program is adopted from the OPTSYS program. Given a 
set of weighting matrices Q and R and the system 
matrices A, B and C, it computes the steady state 
feedback gain F, closed-loop eigenvalues, etc. 

4. Singular Value Analysis Mode : This portion of the 

program allows the designer to analyze various 
designs obtained from the two modes mentioned above 

in term of singular value vs frequency plot. The 

main part of this program is adopted from [Ref. 22]. 

The three modes of operation mentioned above may be used 
in any order to implement specific design objectives. A 
typical design process will involve runs alternating between 
the three modes until a compromise between primary and 
secondary design objectives is achieved. Record of a 

typical design run together with a complete listing of the 

main program and their non-standard subroutine are given in 
Appendix B and C. 

E. DESIGN PHILOSOPHY 

The Linear Quadratic constant state feedback design 
philosophy for the linear time- invariant model is illus- 
trated in Figure 4.2. It assumes that the location of the 
eigenvalues and system time response are the main design 
objectives. These objectives can be achieved using the pole 
placement procedure developed in this thesis. Single or 
multiple reassignment can be made in one run. A number of 
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Figure 4.1 Computer Aided Design Package Organization. 
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designs can then be obtained using different starting 
control weighting matrices, different state weighting 
matrices and different assignment sequences. Physical 
constraints such as control input amplitude, control input 
energy as well as general system properties such as asymp- 
totic behavior are heavily relied upon during this process. 
After the major objectives are satisfactorily achieved, 
secondary design objectives are considered. These include 
feedback gain reduction (by increasing the control weighting 
matrix), robustness (in terms of minimization of system 
sensitivity to modelling errors and/or parameter varia- 
tions), zeros locations, eigenvectors assignments. The 
extra degrees of freedom available in the MIMO state feed- 
back system can often provide a means to improve these 
secondary objectives while only slightly modifying the 
closed-loop pole assignment and thus the time response. 
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Figure 4.2 LQ Eigenvalue Assignment Philosophy. 



V. DESIGN EXAMPLES 



The design procedure described in Chapter 4 is illus- 
trated in this Chapter by two design problems. A typical LQ 
structure and its properties, together with the pole place- 
ment procedure is first presented with a 2x2 model. A prac- 
tical design problem is then presented for the highly 
coupled lateral channels of a CH-47 Helicopter. The 
resulting LQ designs are compared with other multivariable 
state feedback designs [Refs. 22,23]. It is shown that 
the procedure developed here is a viable tool for robust 
constant feedback controller design. 

A. INTRODUCTORY 2X2 PROBLEM 

This problem formulated in reference 23 serves to demon- 
strate how a highly cros s - coupled multivariable control 
problem can be formulated and solved as a linear quadratic 
design problem, using the pole placement procedure. The 
problem provides excellent insight into the structure of the 
multivariable LQ system and its built-in robustness to 
cross - coupled perturbation. 

Figure 5.1 shows a diagram of this basic 2X2 system in 
which the plant is given by 



x(t) = Ax(t) + Bu(t) 



( eqn 5.1) 



y(t) = Cx(t) 



( eqn 5.2) 



where 
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Figure 5.1 A Simple 2X2 MIMO Model. 
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This system has open-loop eigenvalue at -1, -1 and is 
therefore stable. The b |2 term in the control matrix B is 
purposely made large to produce the cross- coupled effect 
from channel two to channel one( see Figure 5.1). The 
design requirement is to select a set of feedback gain F 
such that the closed-loop eigenvalues are at -2, -2. 
Assuming that this is the only requirement, it is not diffi- 
cult to see that a unity feedback law of the form (Figure 
5.2), 



- 
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0 




x > 






u 1 




- 0 


-1 




x , 

L. 






■ _ ^_ l 








L J 




— - 



( eqn 5.3) 

will produce the desired closed-loop system given by, 



- 




















x , 


- 


-2 50 




x , 


+ 


i 


50 




U CI 


X 1 




0 

1 

ro 




x i 




0 


1 




u c2 


J 








_ j 


I 




_J 




— _ 



(eqn 5.4) 

The above design seems to be acceptable as far as eigen- 
values or time response is concerned. It is now shown that 
when robustness of the system is considered, the unity feed- 
back gain controller performs rather poorly. On the other 
hand, design using LQ Pole Placement type of formulation 
will result in robust controllers. To demonstrate the lack 
of robustness of the unity feedback design, the feedback 
gain matrix F as given in equation 5.4 is perturbed slightly 
(by +5%) and the eigenvalues of the resulting closed-loop 
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Figure 5.2 Unity Feedback for the 2X2 Model. 

system matrix are calculated. The absolute and percentage 
errors in the eigenvalue due to the perturbation of the 
closed-loop system are given in Table I . 

It can be seen from the above that the first cut design 
is very susceptible to model and feedback perturbation. A 
5 % changes in the f i( term has resulted in a large shift 
( 160% ) in the closed- loop eigenvalue. Lack of robustness 
in the unity feedback design can also be seen in terms of 
the minimum singular value plot of the return difference 
matrix in the frequency domain. For the unity feedback gain 
system shown in Figure 5.2, the return difference matrix is 
given by, 
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TABLE I 

PERTURBED EIGENVALUE FOR 2X2 MODEL 



Perturbation 
in F (+5%) 


Absolute Changes 
in Eigenvalues 


Percent Changes 
in Eigenvalues 


f , , and f 2 2. 


0.05 , 0.05 


2.5 ,2.5 


f <2 


0.0 , 0.0 


o 

o 

o 

o 


f 21 


0.67 , +3.265 


33.5 , 163.0 


All 


1,255, 5.346 


37.25 , 167.0 



50/s+l 
s+2/ (s + 1) 

(eqn 5.5) 

where G(s) =C(sI-A) ^B is the loop transfer matrix as indi- 
cated in Figure 5.3 

The multivariable Nyquist diagram ( locus of det[I+G(s)] 
for the system is shown in Figure 5.4 If this diagram is 
interpreted as for a single input system, the (-1/2, 00 ) gain 
margin and (±160) phase margin would lead one to believe 
that the design is a good one . This has been shown to be 

not the case, a 5% perturbation in F would cause the system 
to become unstable. The above clearly demonstrate the 

inadequacy of the classical method in evaluating stability 
margin for MIMO system. 

Robustness properties of the unity feedback controller 
will now be analyzed in terms of singular value as discussed 



I + G ( s ) 



s+2/ (s+1) 
0 




Figure 5.3 Loop Transfer Matrix - F = I. 

in Chapter 3. The minimum singular value plot of [I+G(s)], 
is shown in Figure 5.5 as a function of frequency. The lack 
of robustness is clearly indicated by the relative small 
singular value at frequency around 2 rad/s. Using the 
universal phase and gain margin chart developed in 
[Ref. 19], the minimum singular value at this frequency 
corresponds to a gain margin of (0.91, 1.0) and a phase 
margin of (± 4 deg). 

It is now shown that formulation using LQ approach and 
the pole placement procedure developed here will result in 
robust design that meet the time response requirement. 
Furthermore, better insight of the design process can be 
obtained from the procedure to be described here. The first 
step in the LQ design is to determine the asymptotic 
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Figure 5.4 Multivariable Nyquist Plot. 

properties of the system, i.e. movement of closed-loop poles 
as R — * 0 . Using results from Chapter 2, it can be estab- 
lished that as R increase from 0 to oo , both closed-loop 
poles move from infinity on the real axis to the open- loop 
poles location. None of the closed-loop poles stay finite 
as R-> 0 since the dimension of the input control vector is 
equal to the dimension of the state. Assuming that R = I, 
the pole placement is accomplished in two steps. The first 
step is to move the open-loop pole at -1 to -2.0. As the 
system matrix A given is already in Jordan form, no trans- 
formation is required. The pole placement program puts the 
pole at -1.9987 with, 
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Figure 5.5 Singular Value Plots - Unity Feedback. 
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0.00119 0 




0 0 


V 




F ,= 






1 

o 

o 

1 




_ 0. 01983 0 



and 



A 



o ug 



1.9915 0 

0.01983 1.0 



In the second step, the other open- loop pole at -1. is 
moved to -2.0 , using A 0U g as the new plant matrix. The 
resulting Q and F and the augmented plant matrix become, 



-2 . 998 


-149.9 




0.9994 


-49.978 


-149.9 


7499 . 3 


F 2 = 


-0.00841 


0.4211 



The effective Q e and F e required to move both open- 
loop poles at -1.0 to -2. are Q e = Q ( + Q„ and F e = F, + 
Fo as shown in equation 5.6 below. The pole placement 
procedure is completed with the final eigenvalues placed at 
(-1.99255, ±j0. 05628). 



Q e = 



- 


— 






' 


2.99919 


-149.9 




0.9994 


-49.978 


-149.9 


7499.3 


F e = 


0.0114 


0.4211 




_ 




l_ 


— 



(eqn 5.6) 



52 



The singular value plots of the case where R = I 
together with the unity feedback gain ( non-LQ design ) are 
shown in Figure 5.6 The well established fact that LQ 
design possesses ( 1/2 , 00 ) gain margin and (±60deg) phase 
margin can also be readily observed from the same figure as 
the minimum singular values of [I + G] , JL , is greater than 
one for all frequency. Changes in closed-loop eigenvalue 
for a small perturbation in F is again computed as shown in 
Table II . It can be seen that the LQ design is robust with 
the largest percentage change in eigenvalues location of 
only 10%, when compared with the 160% change in the unity 
feedback design. 



TABLE II 

PERTURBED EIGENVALUE FOR LQ DESIGN 



Perturbation 
in F (+5%) 


Absolute Changes 
in Eigenvalues 


Percent Changes 
in Eigenvalues 


f 1 1 and f 2 . 1 


0.2149 , 0.142 


10.7 


, 7.1 


f.2 


0.119 , 0.119 


5.59 


, 5.59 


f 2 1 


0.0592, 0.0578 


2.90 


, 2.89 


All 


0.07 , 0.070 


3 . 52 


, 3.52 



The pole-zero plots of various closed-loop transfer 
functions of the closed-loop transfer matrix for both the 
unity feedback and LQ design are compared in Figure 5.7 and 
5.8 . For the unity feedback design, zero at minus three 
for the input 2 to output 1 channel corresponds to the 
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Figure 5.6 Comparison- Singular Value Plots (2x2 model). 
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minimum singular value frequency ( point A in Figure 5.6). 
It is in fact this zero that causes the the system to be 
sensitive to perturbation. For the LQ design, the built in 
robustness cause the zero at minus three to move toward the 
pole location at -2. 

Improvement in the robustness from LQ design is now 
analyzed in terms of Bode plots. The open-loop gain and 
phase vs frequency plots for a MIMO system can be obtained 
from the open-loop transfer function matrix, G(s). For full 
state feedback, G(s)=F(sI-A) B. (if the open-loop plots 
without feedback is considered, G(s)=C(sI-A) l B ) The 
matrix G is a matrix-valued rational function of s, it 
describes how the system (with or w/o feedback) appears to 
its environment. It is an external description of the 
system and is closely related to the zeros of the system. 
The open-loop Bode plots for the two designs are compared in 
Figures 5.9 to 5.11 As b 2 | =0, there is no coupling from 
channel 1 to channel 2. All channels have a -20 dB/decade 
slope at high frequency which is in agreement with earlier 
observation that each channels has one finite zero. It is 
interesting to note that the unity feedback and LQ design 
result in almost identical gain vs frequency plots for 
direct channels (i.e. channel 1-1 and 2-2). Any classical 
single loop type of analysis will not be able to detect any 
difference between the two designs. On the other hand, 
cross- coupling effect can be readily seen from the gain vs 
frequency plot in channel 2-1 (Figure 5.10). The unity 
feedback design is characterized by the rather large channel 
2-1 gain at low frequency. A very small perturbation in 
channel one's parameter can change the system behavior 
considerably. This has been illustrated earlier by 
perturbing the feedback gain. It can be seen from the 
figure that LQ design reduces the crossfeed gain by large 
amount . 
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Figure 5.7 Closed-Loop Pole-Zero Plots (Unity Gain FB). 
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Figure 5.8 Closed-Loop Pole-Zero Plots (LQ Design). 
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Figure 5.9 Open-Loop Bode Plots Channel 1-1. 
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Figure 5.10 Open-Loop Bode Plots Channel 2-1. 
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Figure 5.11 Open-Loop Bode Plots Channel 2-2. 
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B. A HELICOPTER DESIGN PROBLEM 



The design procedure described in Chapter 4 is further 
illustrated in this section by an actual design problem. A 
controller is designed using the Linear Quadratic Pole 
Placement approach for the lateral dynamic model of a CH-47 
helicopter. The resulting design is then compared to multi- 
variable state feedback controller given in [Ref. 23]. 

The highly coupled two inputs lateral axis model of the 
CH-47 helicopter is used as a full order system. The state 
vector x(t) and control input vector u(t) are given by, 

x, = v = Y-axis velocity (ft/sec) 
x 2 = p = Roll rate (rad/sec) 
x 3 = r •= Yaw rate (rad/sec) 
x + = 4 > = Roll angle (rad) 



u, = = Yaw rate rotor deflection control (inches.) 
u 2 = 8 C = Roll rate rotor deflection control (inches.) 

The state variables and the body axes of the aircraft are 
illustrated in Figure 5.12 The yaw and roll rotor deflec- 
tion control produce changes in the yaw rate, slide slip 
angle, roll rate and bank angle. Assuming full state feed- 
back, the A, B, and C system matrices are given by, 



A = 



-2.27 

.01 

0.04 

0 . 



1.420 

-0.7 

-0.05 

1 . 



-0.15 
-0.07 
-0.05 
0. 11 



31.99 

0.0 

0 . 

0 . 





0.12 


0.95 




1 


0 


0 


0 


B = 


0.04 


-8 . 37 


C = 


0 


1 


0 


0 




.34 


.020 




0 


0 


1 


0 




0.0 


.0 




0 

[_ 


0 


0 


1 
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Figure 5.12 State Variables and Body Axes for CH-47. 
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The open loop eigenvalues of this system are; 

s , = +0.2065 
s ^ = -0.0503 
s 3 = -1.0498 
s 4 = -2.1263 

The open-loop system is not stable and the time response 
is shown in Figure 5.13 for zero input and an initial condi- 
tion of <j>( Q ) = 0.1 rad 

The design requirements are to satisfy specif ication in 
terms of the step input response of the roll attitude 
channel ( 0 / <fi c ). Stability margin requirements are those 
given by standard military specification. In [Ref. 23] 
three designs were obtained to satisfy the desired time 
response performance specification. All three control laws 
are of the form given by, 

u(t) = -Fx(t) + h</> c (t) (eqn 5.7) 

The values of F and h are suminaried in Table III . It was 
shown in [Ref. 23] that two of the designs (design 1 and 2 ) 
were extremely sensitive to model errors and perturbations. 
It is now shown that LQ formulation using the pole placement 
procedure developed here will result in robustness design 
and yet satisfy the conventional time response criteria. 

It is assumed that the closed-loop poles requirement are 
the same as those obtained in the AlphaTech's design 1. 
(-25.12, -12.51, -9.652, -2.125). The first step in the 
design is to establish some asymptotic properties of the 
system. The dimension of the control input (2) is less than 
the dimension of the state(4). As shown in Chapter 2, this 
implies that at least two closed-loop poles approach minus 
infinity in the complex plane when R — » 0. When R -* 0 
(meaning no control input allowed), the closed-loop poles 
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ROLL ANGLE 



OPEN-LOOP RESPONSE 
INITIAL ROLL ANGLE = 0.1 RAD 




Figure 5.13 Open Loop Time Response <£(0) =0.1 rad. 
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TABLE III 

ALPHATECH DESIGNS (1,2, AND 3) 



Designs 


Feedback 


Gains (F) 


h 


one 


-1.72 


-23.5 


70.6 595 


595 




0.024 


-2.71 


0.368 -7.99 


-7.99 


two 


0.198 


154. 


18.3 142.0 


142.0 




-0.01 


-1.592 


-0.189 -1.47 


-1.47 


three 


0 


0 


25.5 0 


0 




0 


-4 


0 -27 


-27 



approach the open-loop poles or their mirror images if they 
are in the right hand plane. (i.e. -0.2065, -0.0503, 

-1.04987 , and -2.12635). 

The next step in the design procedure is to select a 
suitable starting control weighting matrix. Based on guide- 
line given in [Ref. 4] and assuming that u (m0x and u„ mox are 
equal to 1 inch. R = I is selected for the initial design. 
(In the present formulation p is not needed as only Q is 
varied to place the pole ) . 

The sequence of the reassignment is determined next. At 
present there is no known established guideline for 
selecting the prefered sequence of moving the poles. Two 
extreme sequences are considered as follows, 

A. Move the left most open-loop pole to the left most 
desired closed-loop pole etc. The reassignment 
sequence then become: 

1. move pole at -2.1265 to -25.12 

2. move pole at -1.0483 to -12.51 



65 



3. move pole at -0.2065 to -9.652 

4. move pole at -0.0503 to -2.125 

B. Move the right most open- loop pole to the left most 
desired closed-loop pole etc. For the problem given, 
the reassignment sequence are : 

1. move pole at +0.2065 to -25.12 

2. move pole at -0.0503 to -12.51 

3. move pole at -1.0498 to -9.652 

4. move pole at -2.126 to -2.125 

(close, so no move required) 

Using the pole placement program developed in this 
thesis, the corresponding Q and F for each of the reassign- 
ments are obtained and tabulated in Tables IV and V for the 
two cases selected (design LQ-A and LQ-B). It can be seen 
from the table that difference reassignment sequence results 
in different set of Q e and F e . In general, once the 
matrix R and the n closed-loop poles are selected, there are 
n(n+l)/2 extra degrees of freedom available in Q. This also 
verifies the well-known fact that the general solutions of Q 
and R for a given set of closed-loop eigenvalue are non- 
unique. The elements of the Q matrix obtained during each 
reassignment depends on the transformation matrix (M, L or 
U) and hence the eigenvectors that are used to construct 
them. The important of eigenvector type of assignment is 
evident as Q, F and the resulting closed-loop poles depend 
on the eigenvectors used. The designer can shape the 
design by choosing appropriate eigenvectors for the trans- 
formation matrix. The application of these extra degrees of 
freedom will be discussed in the following section. The 
resulting designs (LQ-A and LQ-B) are now compared with the 
first design in [Ref. 23]. 
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TABLE IV 

RESULTS FROM POLE PLACEMENT SEQUENCE (LQ- A) 



Move 


Q and F obtained during each reassignment 




0.10027 


0.96194 


0.11959 


-1.50849 




0.96194 


9 . 22847 


1. 14725 


-14.47182 


Q. 


0.11959 


1.14725 


0.14262 


-1.79908 




-1.50849 


-14.47182 


-1.79908 


22.69429 


F | 


0.00387 


0.03151 


0.00399 


-0.06416 




-0.24260 


-2.83522 


-0.34541 


3 .11642 




6.12853 


-2.99198 


0.51376 


-95 . 14285 




9198 


1.46070 


-0 . 25082 


46.44916 


Q, 


0.51376 


-0.25082 


0.04307 


-7.97589 


» i- 


-95.14284 


46.44916 


-7.97589 : 


.47 7. 05225 


F, 


0.05825 


-0.02844 


0.00490 


-0.90426 




2.27529 


-1. 11079 


0. 19079 


-35 . 32260 




80.37563 


18 .91502 


13 . 64762 


218.27522 




18.91502 


4.45133 


3.21173 


51 . 36732 


Q, 


13 . 64761 


3.21173 


2 .31734 


37.06267 


^ 3 


218.27521 


51.36732 


37 . 06267 


592 .76758 


Ft» 


1.58772 


0.37365 


0.26961 


4.31182 




-8 . 62241 


-2.02917 


-1.46414 


-23.41618 




0.15575 


0.02839 


-2.48989 


0.19453 




0.02839 


0.00517 


-0.45383 


0.03546 


Q, 


-2.48989 


-0.45383 


39.80322 


-3. 10981 




. 19453 


0.03546 


-3 . 10981 


0.24297 


F. 


-0.38007 


-0.06927 


6.07574 


-0.47469 


4 


-0.06411 


-0.01168 


1.02482 


-0.08007 






86.76016 


16.91336 


11.79108 


121 . 81841 


Q e 


16 .91336 


15 . 14567 


3.65433 


83.38010 




11.79107 


3 . 65433 


42 . 30624 


24 .17787 




121.81841 


83 . 38010 


24. 17787 


2092.75684 


F e 


1.26977 


0 . 30745 


6.35424 


2 .86871 




-6.65383 


-5 . 98686 


-0.59394 


-55 . 70242 










2 . 8687 






u(t)= - 


Fx( t) + hi, t , h = 














-55.7024 
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TABLE V 

RESULTS FROM POLE PLACEMENT SEQUENCE (LQ-B ) 



Move 



Q and F obtained during each reassignment 



Q, 


0.00031 

0.05289 

0.00589 

0.04781 


0.05289 

9.06491 

1.00865 

8.19296 


0.00589 

1.00865 

0.11223 

0.91163 


0.04781 
8 .19296 
0.91163 
7.40489 


F, 


0.00017 

-0.01771 


0.02851 

-3.03548 


0.00317 

-0.33776 


0.02576 
-2 . 74347 


Qz 


0.26982 
0.46312 
16 .33746 
10.89106 


0.46312 
0.79492 
28.04222 
18 .69380 


16.33746 
28.04224 
989 . 23975 
659.45728 


10.89106 

18.69380 

659.45752 

439.61426 


F* 


0.44607 

-0.26207 


0.76563 

-0.44982 


27.00978 
-15 .86824 


18.00517 

-10.57837 




0. 16317 
-0.35205 
6.24519 
-9.35123 


-0.35205 
0 . 75956 
-13.47425 
20.17563 


6 . 24519 
-13 . 47424 
239.02527 
-357.90381 


-9 . 35123 
20. 17563 
-357 . 90405 
535 . 90698 


f 3 


0 . 19939 
0.30220 


-0.43018 

-0.65201 


7 .63129 
11.56628 


-11.42668 

-17.31876 




Qe 


0.43330 

0.16396 

22.58853 

1.58764 


0. 16396 
10.61939 
15 .57661 
47.06238 


22.58853 
15 . 57663 
1228 .37695 
302.46509 


1.58764 

47.06238 

302.46509 

982.92603 


F e 


0.64563 

0.02242 


1 0.36396 

1 -4.13731 


34.64423 

-4.63971 


6.60425 

-30.64059 



u(t)= - Fx ( t ) + h <*> c (t ) , h 



6.60425 

-30.64059 



Closed loop time response for a step input in roll atti- 
tude command ( <#> c =0.1 rad) for the three designs are shown 
in Figures 5 . 14 to 5 . 16 . It can be seen from the response 
plots that although all three designs have almost identical 
closed-loop pole location, the step response for various 
states differs. The AlphaTech design's response is charac- 
terized by the large 'overshoot' in the yaw rate step 
response. In the Linear Quadratic Design (LQ-A) , the 'over- 
shoot' occurs with the roll attitude response rather than 
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the yaw rate response; the yaw rate response is well damped. 
The LQ-B design appears to be a better design as coupling 
among modes are small and can not be detected from the time 
response plot. It is noted that the difference in response 
for the same set of closed- loop eigenvalue is due to 
coupling among various modes through their respective eigen- 
vectors. The 'overshoot' in this case is obviously not due 
to complex conjugate poles as all closed-loop poles in the 
three designs are on the real axis. The set of final eigen- 
vector for the three designs are tabulated in Table VI . 
Inter-mode coupling for the yaw rate response in the 
AlphaTech design and roll attitude response in the LQ design 
case A can be readily seen from the table. The issue of 
eigenvector assignment will be discussed in the next 
section. 



TABLE VI 

CLOSED-LOOP EIGENVECTORS 



Eigenvalues /Eigenvector 





-24.79907 


-2.12242 


-9.72403 


-11.9449 


Alpha 

one 


-0.39336 

0.10767 

-0.91306 

-0.00029 


-0.99976 

0.00822 

0.01955 

-0.00488 


-0.48016 

-0.32948 

-0.81181 

0.04306 


-0.41854 

-0.33640 

-0.84283 

0.03593 


1 


-2.12499 


-9.26645 


-12.52066 


-25.20990 


LQ-A 


0.22870 

0.18018 

-0.95603 

-0.03529 


-0.15459 

-0.98218 

-0.01258 

0.10614 


0.01894 

-0.99664 

-0.00396 

0.07964 


-0.12653 

0.99118 

0.00031 

-0.03932 




-25 . 17894 


-9.62708 


-12.51644 


-2.12630 


LQ-B 


-0.12652 

0.99118 

0.00024 

-0.03937 


-0.13209 

-0.98557 

-0.02568 

0.10267 


0.05042 

-0.99375 

0.06075 

0.07886 


-0.99977 

0.00797 

0.01945 

-0.00475 
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CLOSED-LOOP RESPONSE ALPHATEC 1 
STEP INPUT IN ROLL ANGLE 0.1 RAD 




Figure 5.14 Closed-Loop Time Response Plot (AlphaTech 1). 



70 



ROLL ANGLE-RAD 

0.00 0.01 0.02 0.03 0.04 0.05 0.05 0.07 0.08 



CLOSED LOOP RESPONSE LQ-A 
STEP INPUT IN ROLL ANGLE 0.1 RAD 



o is 




Figure 5.15 Closed-Loop Time Response Plot (LQ-A). 
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Figure 5.16 Closed-Loop Time Response Plot (LQ-B) . 
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The difference in time response performance for the 
three designs can also be interpreted in term of the closed- 
loop pole-zero plots in the frequency domain. The pole-zero 
plots for the designs, with <£ c as the input, are shown in 
Figures 5.17 to 5.19 . As all pole locations are the same 
for the three designs, they are indicated only once in each 
diagram . 

In the AlphaTech's design, the zero locations for the 
various channels clearly indicate the weakness in the 
design. Most transmissions zeros are located away from the 
poles locations. The system is therefore strongly coupled 
to its external environment. A good example is the lightly 
"damped" pair of zero at ( - 1 . 05 , ±7 . 3 1 j ) for the v-^ c chan- 
nels. It is in fact these undesired zeros that reduce the 
overall robustness of the system. The mechanism of robust- 
ness improvement in the LQ design can also be seen in the 
pole-zero plots in Figures 5.18 and 5.19 The built-in 
robustness in the LQ design causes the zeros at various 
channels to move to locations where their transmission prop- 
erties can be canceled by the closed-loop poles. An excel- 
lent example is in the v- </> c channel, where the zeros at -2, 
-7, and the mirror image of +24.0 are fairly close to the 
closed-loop pole locations at -2.12, -9.26 and -25.12. In 
the case where the zero from LQ design are not close to the 
closed-loop poles (r- <t> c channel for LQ-A ) the zeros are 
"well damped" and their transmission properties can be 
neglected . 

Robustness properties of the three designs presented 
here can also be analyzed from the open- loop Bode plots. 
The open-loop transfer function gain plots for channel 1-1, 
l - 2, 2-1 and 2-2 for the three designs are shown in Figures 
5.20 to 5.24 Cross coupling problem for the AlphaTech 
design is clearly indicated by the relatively high gain of 
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F igure 5 . 17 



Closed-Loop Pole-Zero Plots (AlphaTech 1). 
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Figure 5.18 Closed-Loop Pole-Zero Plots (LQ-A) . 
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Figure 5.19 Closed-Loop Pole-Zero Plots (LQ-B). 
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the channel from input 2 to input 1 (2-1). Both LQ designs 
presented here reduced the gain in this channel by more than 
40 dB in the frequencies of interest. The bandwidth is 
reduced from above 100 rad/s to about 8 rad/s. Unlike the 
simple (2x2) system where there was little difference 
between the direct channels of the two designs, gain adjust- 
ment is observed in all channels. For example, LQ designs 
reduce the gain in channel 1-1 but increase the gain in 
channel 2-2. There is also a slight increase in the 
coupling channel 1-2. The overall effect is that of gain 
balancing, gains in channels that are affected by cross- 
coupling perturbation are lowered together with some adjust- 
ments in other channels. Different reassignment sequence 
results in different adjustments. The designer has to 
choose a set of gain curve depending on the particular 
requirement. The relationship between the open- loop Bode 
plots and the zeros of the system is also evident from these 
diagrams. The low frequency resonance (or 'peak') in the 
gain vs frequency plot for the AlphaTech design correspond 
to the undesirable zero mentioned earlier. In the LQ 
designs, these low frequency resonances are absent because 
of the more desirable zero locations. It is also inter- 
esting to note that the LQ-B gain curves fall nicely in 
between the other two designs. This fact, together with its 
better overall low frequency gain characteristic , may 
account for the better time response behavior of the LQ-B 
design . 

As a final comparison, the closed-loop frequency 
response plots of the various channel for the three designs 
are shown in Figures 5.25 to 5.32 . At high frequency, all 
the gain plots approach the -20dB/decade slope as all eigen- 
values are on the real axis. The AlphaTech design has a 
near Odb gain for frequency up to about 2 rad/s. It appears 
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Figure 5.20 Bode Plots Comparison- Input 1-1. 
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Figure 5.21 Bode Plots Comparison - Input 1-2. 
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Figure 5.22 Bode Plots Comparison - Input 2-1. 
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Figure 5.23 Bode Plots Comparison - Input 2-2. 
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to be an acceptable design but was shown to have undesired 
transmission zeros in other channels; this can also be seen 
from the closed- loop plot for the r-«£ c channel in Figure 
5.31 . Design LQ-A, obtained from the reassignment sequence 
A, is characterized by very low DC gain (-9dB) and 'peak' in 
frequencies near 10 rad/s. These have been shown, both in 
the time response and pole-zero discussions that it will 
have undesired effect on the pilot’s control. Design LQ-B, 
obtained from the second reassignment sequence appears to be 
the best compromise. All channels (Figures 5.30 to 5.32) 
have flat low frequency characteristic and coupling between 
modes are almost absent. 

C. DISCUSSION AND CONCLUDING REMARKS 

A new computer aided design procedure for the multivari- 
able linear time- invariant system using Linear Quadratic 
Pole Placement formulation is presented. The two design 
examples presented above have served to demonstrate the 
complexity involved in a multivariable design. The main 
problem lies in the fact that the solution of a MIMO problem 
is in general non-unique. This was shown in the helicopter 
problem where different approaches result in different 
designs, although the closed-loop eigenvalues for all 
designs were the same. It was also shown that the extra 
degrees of freedom in MIMO system design can be accounted 
for by analyzing the singular value plots, transmission zero 
movements and closed-loop eigenvectors of the designs. It 
becomes apparent that the success of any multivariable 
design methodology hinges on how to make use of these extra 
degrees of freedom available in MIMO system. It was also 
shown in the above design examples that the Linear Quadratic 
Pole Placement procedures possess such quality. Some unique 
properties of the method presented here are as follows; 
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Figure 5.24 Closed-Loop Bode Plot- AlphaTech 1-1. 
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CLOSED-LOOP BODE PLOT ( r -</>.) 



INPUT # = 1 

OUTPUT # = 3 
TF GAIN = 2.021*10 




Figure 5.25 Closed-Loop Bode Plot- AlphaTech 1-3. 
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Figure 5.26 Closed-Loop Bode Plot- AlphaTech 1-4 



CLOBED-LOOP BODE PLOT ( V - 9p 



INPUT # = 1 

OUTPUT # = 1 

TF GAIN =-5.2 57 • 1 C 




Figure 5.27 Closed-Loop Bode Plot- LQ-A 1-1. 
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MAGNITUDE - DECIBELS 



CLOSED-LOOP BODE PLOT ( r ~ 4 > c ) 



INPUT # = 1 

OUTPUT # - 3 
TF GAIN =- 1 . 3 8 6 • 1 ( ' 




Figure 5.28 Closed-Loop Bode Plot- LQ-A 1-3. 
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MAGNITUDE - DECIBELS 



INPUT # = 1 

OUTPUT § = 4 
TF GAIN = 4. 664*10 



CLOSED-LOOP BODE ( <t> - ) 




Figure 5.29 Closed-Loop Bode Plot- LQ-A 1-4. 
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Figure 5.30 Closed-Loop Bode Plot- LQ-B 1-1. 
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MAGNITUDE - DECIBELS 




Figure 5.31 Closed-Loop Bode Plot- LQ-B 1-3. 
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MAGNITUDE - DECIBELS 



CLOSED-LOOP BODE PLOT ( <t> - 4> c ) 



INPUT # = I 
OUTPUT f = 4 

TF GAIN = 2 560*10 




Figure 5.32 Closed-Loop Bode Plot- LQ-B 1-4. 
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1. It guarantees stability automatically. 

2. It has partly overcome the main difficulty in 
applying optimal control theory to practical prob- 
lems; i.e. that of selecting suitable performance 
indices (Q and R matrices). The designer can now 
choose to compute Q (if he knows somethings about R) 
or vice versa or compute both Q and R if both of them 
are unknown. 

3. It has built-in robustness to model's error and 
perturbation . 

4. Reduced-order type of problem formulation is possible 
as the procedure can reassign poles one at a time 
without affecting others. 

On the other hand, the design procedure developed here 
is not without some shortcomings. Three main problem areas 
are described below, one of which can be overcome by addi- 
tional programming efforts while the others are still active 
research areas pursued by many researchers. 

1. Eigenvectors Assignments: It was shown in the design 

examples presented above that while the overall speed 
of response of the closed-loop system is determined 
by its eigenvalues, the 'shape' of the transient 
response depends on the closed- loop eigenvectors. 
The problem of eigenvectors assignments was first 
formulated in [Ref. 16]. Since then many eigenvec- 
tors assignment algorithms have been developed 
[Refs. 17,24] for use in multivariable design. In 
principle, the eigenvectors assignment routine can be 
incorporated into the coordinate transformation 
portion of the OPTPP program (as indicated in Figure 
4.1 in Chapter 4). In essence, a new similarity 
transformation different from M in Equation 4.6 is 
computed once the eigenvalues and eigenvectors are 
specified. As most eigenvector algorithms available 
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at present are iterative in nature and their inclu- 
sion requires major programing efforts, it is recom- 
manded for future work. 

Perturbation or Model's Error Structure: The guaran- 
teed robustness obtained from the LQ formulation 
given by equations 3.4 and 3.5 ensures that the 
perturbation or model's error (L(jw)) is sufficiently 
small so that the closed-loop system remains stable. 
However, the above only applies to simple model's 
error structure where both L(jw) and R are diagonal 
matrices. There may be cases when the above equa- 
tions do not hold and hence the design becomes very 
conservative. An example of this nature is shown in 
[Ref. 19]. An ad hoc solution is to use non-diagonal 
control weighting matrix as mentioned in Chapter 4. 
Two designs for the helicopter problem are obtained 
using non-diagonal R. Their results in terms of 
singular value plots, Bode plots etc are compared 
with other designs in Appendix A 

Reassignment Sequence: As shown in the last section, 
different reassignment sequences result in different 
designs. At present there are no known methods for 
determining the best reassignment sequences. 
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VI. CONCLUSIONS 



It was demonstrated that the general pole assignment 
problem in multivariable state feedback control system 
design can be formulated using the Linear Quadratic Control 
approach. This method of formulation is effective for two 
main reasons; First, the extra degrees of freedom available 
in a multivariable system structure is utilized to produce 
designs that are robust to perturbations in the system and 
gain matrices. Secondly, the classical difficulty of 
selecting suitable performance index in optimal control 
formulation was partly overcome, as designers now have the 
flexibility of specifying only Q (the state weighting 
matrix) or R (the control weighting matrix) or both Q and R 
as the design parameters to be varied. In other words, 
knowledge of the performance index which ideally should come 
from physical argument is used to the best of designer's 
advantage. In addition, the structure of the present formu- 
lation is such that eigenvector assignment can be further 
incorporated into the procedure. The above properties, when 
combined with the reduced-order formulation capability, have 
been shown to be very versatile and have important impact on 
the performance of the resulting design. 

The optimal pole placement (OPTPP) program developed 
here is combined with other well established routines to 
form a computer aided design and synthesis package. 
Together with the design procedure and philosophy presented 
here, it provides the control system designer an excellent 
and viable tool to solve complex multivariable problems. 

The procedure was applied to practical test examples, 
and numerical results were presented and discussed. Results 
indicated that all controllers obtained from the formulation 
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given here were stable and robust. Introducing perturbation 
in the system matrices leaded only to small errors in the 
assigned poles. The main shortcoming of the design proce- 
dure is the ad hoc nature in which the poles are reassigned. 
More research is required to develop a systematic way of 
assigning poles and its eigenvectors thus allowing the 
designer to optiinaly shape the response of the system. 
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APPENDIX A 



NON-DIAGONAL R DESIGNS 

The two design examples presented in Chapter 5 were 
based on diagonal control weighting matrix (i.e. R= I). It 
is now illustrated that design using non-diagonal R will 
provide yet another degrees of freedom available to the 
designer. This type of formulation is especially useful 
when the model's or error structure is known. For multipli- 
cative type of perturbation, the effect of control weighting 
matrix on the system stability has been explored in 
[Ref. 12]. In essence, the selection of the R matrix deter- 
mines the coordinate frame in which the sensitivity assess- 
ment is to be made. This can be readily seen from equation 
3.1, the general form in which R is non-diagonal is given 
below; 

1/2 —1 -1/2 

« [R L (s)R - I] < 2. [I + G(s)] ( eqn A.l) 

In general, the sensitivity of the system to perturba- 
tion in a particular L(s)'s direction can be reduced by 
making a [R ^ L (s)R - I] small. This can be done 
simply by choosing an appropriate R. The main problem with 
this kind of approach is that L(s) must be known precisly 
for all frequency for the computation of <T [r'^L* 1 (s)R - 
I]. In addition, the worst case direction of L(s) must be 
known otherwise the resulting design may be too 
conservative . 

The helicopter problem presented in Chapter 5 is now 
analyzed using the non-diagonal control weighting matrix. 
Two designs are obtained as follows; 
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Design one (LQ-C): 



1.0 


0.0 


0.0 


2.0 



Design two (LQ-D): 



R 



1.0 

0.5 



0.5 

1.0 



Design LQ-C is obtained using an assignment sequence 
similar to LQ-A with no off-diagonal element in R. Design 
LQ-D is obtained using placement sequence for design LQ-B 
but has off-diagonal element in R. Q and F obtained during 
each move and the final Q e and F e are tabulated in Tables 
VII and VIII. Singular value plots and the open-loop Bode 
Plots are compared with other designs in Figures A.l to A. 5 
Both design has singular value greater than one for all 
frequency. It is interesting to note that similar reassign- 
ment sequence produce similar singular value plots. This 
can be readily seen by comparing LQ-A with LQ-C and LQ-B 
with LQ-D in the figures. The effect of using non-diagonal 
R merely changes the shape of the singular value plots. For 
the cases presented here, both non-diagonal R designs are 
slightly less conservative (having lower singular value). 

As in the singular value plots, the shape of the open- 
loop Bode plots (with feedback) are closely related to the 
reassignment sequence .( compare Bode plots of LQ-A with LQ-C 
and LQ-B with LQ-D in Figures A. 2 to A. 5 . Results from 
pole-zero maps also indicate similar trends. 

In summary, it is demonstrated that designs using non- 
diagonal control weighting matrix (including possibly off- 
diagonal elements) provide yet another means of 'fine 
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tuning' the design. This capability, such as using the off- 
diagonal elements directly in design, is unique to the 
present formulation. Robustness between the upper and lower 
crossfeed (of L(s)) can be controlled by adjusting the rela- 
tive weighting of the upper and lower (or off-diagonal ) 
elements of the control weighting matrix. It must be empha- 
sized that this kind of fine tuning is only possible for a 
class of rather well-defined structure of L. In practices, 
other constraints, such as the energy of the control input, 
the conditioning of R etc must be considered over the range 
of frequencies. The key issue of model's error structure 
and how it can be used in multivariable control system 
design is currently being investigated by many researchers. 
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TABLE VII 

RESULTS FROM POLE PLACEMENT SEQUENCE (LQ-C) 



Move 



Q and F obtained during each reassignment 



Qi 


0.00062 

0.10604 

0.01180 

0.09584 


0.10604 
18 . 17288 
2.02209 
16.42485 


0.01180 

2.02209 

0.22500 

1.82758 


0.09584 

16.42484 

1.82758 

14.84496 


F, 


0.00033 

-0.01773 


0.05708 

-3.03892 


0.00635 

-0.33814 


0.05159 

-2.74657 


Qi 


4.17258 
-3.78394 
0.13594 
■108 . 79944 


-3 . 78394 
3.43150 
-0 . 12328 
98 .66574 


0. 13594 
-0.12328 
0.00443 
-3 . 54455 


-108.79945 
98 . 66574 
-3 . 54455 
2836.93335 


F z | 


0.02701 

1.21663 


-0.02449 

-1.10329 


0.00090 

0.03967 


-0.70429 

-31.72316 




154.23135 
36.25874 
26.52100 
418 .56201 


36 . 25874 
8 . 52419 
6.23491 
98.40114 


26 .52098 
6.23491 
4.56044 
71.97421 


418 .56201 
98.40111 
71.97421 
1135.91846 


f 3 


2.80842 
-7 . 58133 


0.66025 

-1.78234 


I 0.48291 I 

1 -1.30355 


I 7.62180 

1 -20.57501 


<4 


0.44288 

0.09189 

-4.22665 

0.82784 


0.09189 
0.01906 
-0.87692 
0 . 17176 


-4.22665 
-0.87692 
40 .33727 
-7 . 90052 


0.82784 
0. 17176 
-7 . 90052 
1.54741 




-0.63449 

-0.09950 


I -0.13165 
1 -0.02064 


6.05549 

0.94961 


-1.18605 

-0.18599 


F e 


2.20127 

-6.48193 


I 0.56119 

1 -5.94519 


6 . 54565 
-0.65241 


5 . 78305 
-55.23073 





158.84741 


32.67273 


22.44206 


310.68604 


Q . 


32.67273 


30.14761 


7 .25680 


213.66344 


^ e 


22.44208 


7 .25680 


45.12712 


62 . 35670 




310.68604 


213.66348 


62 . 35670 


3989 . 24390 



u(t)= - F x ( t ) ♦ h 4> c ( t ) , h 



5 . 7839 
-55.23 
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TADLE VIII 

RESULTS FROM POLE PLACEMENT SEQUENCE (LQ-D) 



Move 



Q and F obtained during each reassignment 



Q. 


0.00023 
0.03940 
0 . 00438 
0.03561 


0.03940 
6.75195 
0. 75 L29 
6 . 10249 


0.00438 

0.75129 

0.08360 

0.67902 


0.03561 
6 .10249 
0.67902 
5 .51550 


F, 


0.00895 

-0.01765 


1 .53374 
-3.02506 


0 . 17065 
-0.33659 


1 . 38619 
-2.73405 


Qz 


0.24621 
0 . 11866 
14 . 98355 
9 .73272 


0. 11866 
0.05719 
7.22105 
4 .69050 


14 . 98355 
7.22105 
911.83813 
592.29370 


9 .7 3271 
4 . 69050 
592.29370 
384 .73022 


Fz 


0.56910 

-0.32115 


1 0.27424 

-0. 15476 


34 .63293 
- 19 . 54404 


22.49620 
- 12 .69504 


Q3 


0. 17288 
-0.43406 
7 . 14036 
-7.87012 


-0.43406 
1 . 08987 
-17 . 92834 
19 .76067 


7 . 14036 
-17.92834 
294 . 92139 
-325.06299 


-7.87012 

19.76067 

-325.06299 

358.28516 


F J 


0.05761 

0.34047 


-0 . 14465 
-0.85487 


2 .37951 
14 . 06264 


1 -2.62273 

| -15.49990 



F e 


0.63566 

0.00167 


1 . 66333 
-4 . 03469 


37 . 18307 
-5 .81798 


21.25964 

-30.92899 




Qe 


0.41932 
-0.27600 
22 . 12828 
1.89821 


-0.27600 

7.89901 

-9.95601 

30.55365 


22 . 12828 
-9.95601 
1206 . 84302 
267 . 90967 


1 .89820 
30 .55365 
267 . 90967 
748 .53076 



u(t)= - Fx(t) + h * c (t) , h 



21.159 

-30.9289 
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SINGULAR VALUE PLOT 



o 




Figure A.l Singular Value Plots - Comparison. 
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MAGNITUDE - DECIBELS 

40.0 -20.0 0.0 20.0 40.0 GO.O 00.0 



OPEN LOOP GAIN 1 - 1 




Figure A. 2 Bode Plots Comparison- Input 1-1. 
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MAGNITUDE - DECIBELS 

60.0 -40.0 -20.0 0.0 20.0 40.0 60.0 



OPEN LOOP GAIN 





10 10 10 
FREQUENCY (RAD/S) 



Figure A. 3 Bode Plots Comparison - Input 1-2. 
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MAGNITUDE - DECIBELS 

40 00 120 100 200 



OPEN LOOP GAIN 2-1 




Figure A. 4 Bode Plots Comparison - Input 2-1. 
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MAGNITUDE - DECIBELS 



OPEN LOOP GAIN 2-2 




Figure A. 5 Bode Plots Comparison - Input 2-2. 
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APPENDIX B : EXAMPLE DESIGN RUN OUTPUTS 
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APPENDIX C : COMPUTER PROGRAM LISTINGS 



ooooooo 

~ 4 (%j m-rmsOh- 



UUUOOQO 

ooooooo 

Z 2 Z 2 ZZZ 

UOUQOOO 

UOUUUUO 



oooooooo ooooooo ooo ooooo 

o-Hcgi^-si-^-oi^* coO'O^r^rO'-t LAsOr^ coo'o^cvl 

i\|iNji'jLNi' 4 L\| , 'JC\giN 4 c\imcncMCi|rncMrnrT) nn'i'T'T'T 

ooo ooooo ooooooo OOO OOOOO 

OOOOOOOO OOOOOOOOOO ooooo 

2 ZZZZ 2 ZZ 22 Z 2 Z 2 ZZ 22 2 ZZZ 2 

□uuGooauuDuoaoouoa qquou 

ooo ooooo ooooooo ooo ooooo 



ft ft ft 



ft 

ft 

ft 

ft 

z 

< 

CC 

o 

o 

oC 

CL 



UJ 

2 - 

LU 

o 

< 



LU 

-J 

o 

CL 



< 

X 



CL 

o 



CL 

CL 



ft 

ft 

ft 



ft ft ft ft ft ■ 


ft ft ft ft ft ft 


*###***# 


ft ft ft ft 


ft ft ft ft ft ft ft ft ft ft ft 


ft ft ft ft ft 




X 




-* 




tt 










0 


►~<LU 


a 


tt 


• • 


ft 










• UJ 


OCcC 


LU 


# 


0 


ft 




UJ 






CL h- h- 


H-< 


CC 


* 


0 


tt 




h* 






XOUJZUJG 


< 


►— « 


# 


a 


ft 




< 






—oxujh-CL 


zz 


O 


ft 


CO 


ft 




h- 






cLjoroz 


LU 


O 


ft 


X 


ft 


cC 


O 




z. 


h 1 rrcLO 


OX 


UJ 


ft 


> 


ft 


ox 


z 




a 


< 0002:0 


Xh- 


X 


ft 


0 


ft 


z 








XUJX.-.0 


*— • 




ft 




ft 


O ^ + 


X 




Zh 


OhOUZ 


h- LL 


X 


ft 


ox 


ft 


—• O 






o< 


ooxo uj 


xa 




ft 


z:< 


ft 


h- <c>- 


H 


• » 


►— 0 


ZJO<OX 


0 




ft 


< CC 


ft 


CL OX 


ZD 


LU 


1 - h-a 


-•O'— | — h* 


— X 


0 


ft 


0 


ft 


*-« Q 1 


CL 


h- 


ZD <0 


1— OUJ 


U J h- 


UJ 


ft 


JO 


ft 


CC UJ 11 


h- 


< 


CL O 


XOOJQIL! 


'SO 


1 — 


ft 


<X 


ft 


0 LUO 


ZD 


h- 


h- OUJ 


Ouj<olux 


CO 


< 


ft 


LUCL 


ft 


0 X 


OO 


in 


3 00 -J 


Md o j < 


0 


a: 


ft 


OC 


ft 


UJ CC 


oz 


2 


0X0 


UJMX 0 


OX 


O 


ft 


LL 


ft 


Q..LUQ-. 


a: — 




"Sh^OUJCL 


soox <0 


uC 


a. 


ft 


OO 


ft 


Oh- UJ h- 


X 


h-h-»-H O 


UJ OUh 


h- 


oc 


ft 


LU 


ft 


x<<xi— ZX 




JGhQO 


LUOIU O 


XX 


0 


ft 


QC|- 


ft 


UJ h- + <0~ 


_j 


CL CL < CL UJ 


1— O X O — ! 


Oh 


0 


ft 


<0 


* 


HZO h-OOC 


0 


zzu a: 


<uj *— * O • Cl 


0 X 


21 


ft 


JH 


ft 


XUJ XOZh 


a: 


*— • mO 0 ♦ 


HIQh 


h* 


►—i 


ft 


00 


ft 


co>xxz: 


< — 


h-^ 


'-'-jyjo 


Oh XI— UJUJ 


Ul< 




ft 


UJ 


ft 


—0 1 


xxo^o 


o:uj 


UUJZG 


XX 


UJ 


ft 


Q< 


ft 


aox 11 x 


zoz 


XOmQ 


luz: 0^0 




CO 


ft 




ft 


2 X 


UJh-UJOLU 


•— *uj 0 


X LU U i CL 1 — < 


O 




UJ -fr 


OO 


ft 


<0 «* OJhGJZJXCLJUJUJ 


h>IZZ> 


luz: 


z: 


CC ft 


Zh 


ft 


•-» 2_l— < A. CD 


^CQ^hOOI 


hHOG 


h- ►— • 


< 


0 ft 


< 


ft 


UJ CC<hhOXO(L<CL h- 



LUO OC QC • 
h- O X <0 
•J 2 - 1 Z -JO 
diuIOuoGUJ 

2:0 >OZ 

ocoa 02 : 1 — 

OO UJ H-.U0 
Q1ZZQ.L03 
O CL 0^)0 CO 



Dh 
CL X 
SO 

o~ 

Oil! 

S 

Q 



Ul 

2 : 
2: 
o 
a 



CL 

< 

o 



O 

LL 

O 



Sf 

-* 

ft 

# 

* 

# 

ft 



o 

ooz 

x~ 

os 

mQ 

O-J 



ujO 

XX 



■tf- 

-ft 

tf 

tt 

ft 

ft 

# 



h- *-»0C 0 


oa: 


0 






*— 1 


* 


-JO UJ • CrC 


Ul |— 


0 


OLA 


ft 


OUJ 


ft 


QOUJUJ XX 


•— . z 


< 


— '00 


ft 


X 


* 


LUc^OXl— »-«x 


xo 




-J JN 


ft 


Oh 


ft 


Zh O OCLU 


— 'O 


X 


aj^ 


ft 


h~ 


ft 


OZUJ XI— h- 


0 


0 


O 


ft 


UJ • 


ft 


h-oxoo<o 


00 




CL h 


ft 


OO 


ft 


OUhZ 2! >- 


XX 


UJ 


UJCL 


ft 


zz 


ft 


LU >~+Z O 


< 


X 


CCUJ 


ft 


^0 


ft 


Oo alOZ 


Ui 


> 


0 


ft 


Oh 


ft 


< OmmIIJ 


CO UJ 


h- 


9* 


ft 


Ol- 




0 rz 0 O O < X 


h- 




0 


ft 


OO 


ft 


*hhZ CCOH 


X< 


CC 


X • 


ft 


oj^ 


ft 


i — O OUJ 


<h- 


0 


UJ **"4 


ft 


O 


ft 


ZQ.mUJ>^UJ 


OO 


h- 




ft 


h-X 


* 


< Oh 1 O Z 




O 


z 


ft 


0 




UL <IUQ<1^ 


XX 


X 


xo 


ft 


ax 


1C 


UOUdUJLUi. 


<x 


J> 


<►-* 


ft 


—X 


ft 


GZQ MQd 


CC h~ 


X 


CO 


ft 


u 1 


ft 


a:— foujxuiuj 


O • 


UJ 


q: 


ft 


CL 


ft 


U-il (X^LLih- 


OIZ 


O 


-CLU 


ft 


U)2I 


ft 


oooxui 


OCh- JC 


»— » 


a> 


ft 


* — * C — ) 


ft 


o<ujZg q 


CLOO 


LU 


X 


ft 


xo 


ft 


H-. 0 XLU 

iiiOd xa 


COX 

lu 


X 


0 


ft 

ft 


h- 


ft 

Sc 


h-oa-o<i-i— 


XQCX 

hQG 


X 

h- 


BY 


J* ft ft ft ft ft ft 

ft 

ft 



STS! OH-OOX^ah-Z: UJ 

<axii^zoaixcL< axu. 
z~ft 2 :^ 1 - zaoho 

hOUJ^ XXC*C ZG 

X X GUZuO^Jllh 

CD a: I! I-XX < <Xh- I Oa 
U J >- ►-»*-» *03X0*— »XO < 

OO xxa>o^ai— oujh- a. 

UJ*-H ••OHhd LLKUiX'—' OX 
OO <<LUUJOUJOUJO<X 
CQZX^ZZ^li. LL ►— « J2 _J CL UJ 

<0 O CD X UJ O I 

^U^«JGLHUlll^J _J CL 

oc co o h-<»— aa<z 

<CC-fQh-0Ch-<C0<LUCXUJUJO 
>UJ LUXI— ZHQh-H h-XXO 
QXUJ<Z<uOLU;/l<Z^ 
UZiHLJO J IU h-OOUJUJ 
<Z)< CLOQ. XXXOOUiXX 
LU h- OO 01— H-OLL 

CXl'Fl II OUILUUJU-UJLLUJUI ft ft < 

UJ dill X XXUJ | | LL CL) 

I 



o 

ft 

X 

1 C 

CO 



o 

ft 



000 
>-00 
ox I 



<COO X 



h- h- X 

h-OZD 

Z 2 I I O 
I I I OO X 

< < O CO 
220 LUZOLLZ 
Oarsiai'— llco< 






Sf- 






o 000 00 00 00 000 00 00 0000 000 000 00 000000 ooooooo 00 000 



116 



o oooooooooooooooooo ooooo o oooooo 

rn mvOHcoo o ogm>rtn'OM 

>j* >TNr'Tu\umtutu\a\u\u\u\u^vU'4j^vu r*- i^r^p*HP*r» 

o uoooaouuaouoouaQoo ooooo o oooooo 

o oooooooooooooooooo ooooo o oooooo 

2 Z2222Z222ZZ222222Z ZZZZZ Z ZZZ2ZZ 

a uuauauouooouaououo aouou a aoaaoa 

o oooooooooooooooooo ooooo o oooooo 



o 

ON 

I'* 

o 

o 

z 

o 

o 



ftftftftftftftftftftftftftftftftftftftftftft 



ftftftftftftftftftftftftftftft 









1 


— 


ft 




ft 




1 1 








LU 


1 




ft 




ft 




1 ^ 1 








Q 


1 


-% -*• X 


ft 




ft 




1 O 1 






1 


0 


1 


to — * 0 


li- 


•» — % 


ft 




1 O | 






1 


Z 


1 


— 0 — 


ft 


— O 


ft 




1 O | 










1 


Z LL > 


li- 




ft 




1 CD | 








00 


1 


O — Z 


ft 


0 '— <* •» 


ft 




1 1 








►— 1 


1 


*— 1 • — % > •— « 


ft 


0 0-* 


ft 




1 Z 1 










1 


1 - OO X Z + 


li- 


*0 0 00 


ft 




1 O 1 








> 


1 


O H ft >«H 


ft 


— *in ^rH 


ft 




1 M 1 








-J 


1 


Z 3 O + — 


ft 


OUJ-.— CC~ 


ft 




1 H 1 








< 


1 


X Q- >-1 


ft 


in^roo a.o 


ft 




1 < 1 








z 


1 


LL Z + — LULU 


li- 


— * *0 cC 


ft 




1 2: | 








< 


1 


^ — LU 33 


ft 


3 -^LOCL U 1 * 


ft 




1 DC 1 










1 


DC Jh ugjj 


li- 


ro-o ho 


ft X 




1 O 1 




2: 




UJ 


1 


LU • < — 3 -J<f < 


ft 


OOO h *z 


ft <r 




1 X | 




LU 




X 


1 


LL O — >> O <t>> 


ft 


< 0 >hoO — 0 


ft CC 




1 oO | 


— * 


1 — 




-J 


1 


00 -*z niai<> 


ft 


X H LU * OO 


ft 0 




1 x 1 




OO 




< 


1 


Z — -*OftCOO> CCCC 


li- 


* hq C 


ft 0 




1 < 1 


00 


> 




> 


1 


<— ll 11 o^iizjj oc<< 


ft 




ft DC 




1 CC 1 


t—t 


00 






1 


a:-* ft ll h<<io:<jj 


ft 




ft CL 




l 1- 1 


00 






a: 


1 


H ftoOX +00>><JOO 


ft 


0 * -mo xo 


ft 




1 | 


u.t 


LU 




< 


1 


O H + 0300 


ft 


LO-— — • ^LT\ ^ 


ft LU 




1 X | 


X 


X 




0 


1 


H *^ 00 + ^OoCCCiOO^Z 


ft 


^oaX'- 00 * 


ft X 




1 H 1 


H 


H 




X 


1 


Zft LL— a. mwZ«(DZ^ m 


ft 


3ooxx * *0 


ft H 




1 < 1 








0 


1 


<LL H>-«^ OOOZ»"*oOOO 


ft 


3; *0 


ft 




1 x | 


X 


x 






1 


— 1— ft >3 — lloXXhhOO 


li- 


jO'-OZOOv 


ft z 




1 * —i 1 


X 


O 




>— ♦ 


1 


CL> ZQ LLOLUOO00 hh 


ft 


< 0 <>H> 0 ^ f H 


ft H -1 




1 a 1 


H 






00 


1 


ZO^ LL O </)ZZ t -33 


-* 


UJ— s 0000 in — CL 


ft 




1 CC | 




UJ 






1 


ft •-« *0 LU — H 3 Q_Cl 


ft 


oco * ^ *—000 


ft 0 




1 cj 1 


Z 


X 




LU 


1 


+ +Q Lu 3 ft 000030 .HH 


ft 


*LL-^^— <N<X 


ft LU 




1 0 1 


►—1 


O 




X 


1 


LL + ZUOJ Q.Z 33 


ft 


^rQOOorvj>a: 


ft 00 


— « 


1 0 1 




< 




H 


1 


•— M“H JJ<JhhZHQJ 


ft 


OX— <000.00 * 


ft O 


C\J 


t 1 


X 


> 






1 


+ *-< LLJOOOh 


ft 


o< *lt\lo 0 *— 


ft 


N 


1 x | 




z 




>- 


1 


T>CL Cl lulliuj 


li- 


► *0— *-^o 


ft 00 




1 X 1 




LU 




f£) 


1 


*-* > QC h-HUJ>>> 


ft 


O— 4 CNJX— O-l 


ft aC 


1 


1 H 1 




O 






1 


— XXXH QMODO>hmh 


li- 


oO — Z>-iOo ^ 


ft UJ 


w 


1 1 


< 


•— « 




Q 


1 


aouooct:<jzoa^hhh 


ft 


— o^ZXO—O 


ft H 


ft 


1 z 1 




LU 




LU 


1 


LL LU<_ 1 D M H<<< 


li- 


0 ^ so>inMj 


ft UJ 


23 ft 


1 1 


X 






00 


1 


OLUUJLUa:o 3 L 000 UiIU<OOO 


ft 


*0 *^oO^-X — 


ft X 


a tcccc 


1 1 


0 


O 




3 


1 


333 X 30 Z »Ommm 


li- 


— ^<^»o-^oO *0^*0. 


ft c 




1 O 1 




LU 






1 


LU J J Jh O Z^t-MHH JJJ 


ft 


MO— O *^C\JOOOO 


ft OC 


XXX 


1 x 1 


DC 


H 




LL! 


1 


D<«tzz^oooohh-iaaa 


ft 


| oUH^OQ. 


ft <<COOLLOOO 


1 00 1 


O 


X 




DC 


1 


> J>>>m»-hLO DC *— » • CL *—• •— 1 >-h 


ft 


O Hi. — OOJ — CC 


ft o. 




1 3 1 


H 


cl 




<T 


1 


o<r 00 ZMQQhhhH 


ft 


*om >oin *0 * 


ft xxxxxxx 


1 1 


O 


i 






1 


z>a:a:a:x qdilqqjjjj 


ft 


X— • *oOLO — 


ft >< LULU LULU UJU.I LU 


1 X 1 


Ul 


n 




00 


1 


LU <«OhZr^<ODOD 


ft 


| — ^wQoi^O 


ft LJLJ_I 1 — l 


__j » 1 1 


1 DC 1 


> 


0 




O 


1 


DoOJJ l/) 0*— *h* ZZZZ 


li- 


<codo — oxo — — 1 


ft — IQ.CLO 


Q_ CLCL CL 


1 < 1 


z 






z 


1 


o<ddo>q:oxo:zx 


ft 


— *— IOXOLOO * 


ft Q.XXXXXXX 


1 1 


X 


0 




►— < 


1 


iujoooj^iu<a*H<zxzx 


ft 


n**^ *mz< w xo 


ft 20000000 


1 00 | 


0 


< 




3 


1 


onzzzzu,^ 5 'oozz*H<M< 


li- 


ft 00 — 0>MXo 


ft 00000000 


1 O 1 




LU 




0 


1 


llom^mo triz 


ft 


JJHOh(/)H>w 


ft 0 




1 z 1 


LU 


DC 




— J 


1 


I^t/)oOoO 1 1 1 1 1 1 


ft 


<1 *— *— <00 *U.OOJC 


ft III 


1 1 1 1 


1 ►“» 1 








X 


1 


1 -h 1 1 1 1 


li- 


LL 0 O UJ X -CL 


ft UJ 




1 2 1 


1 


1 




a 


1 


00 1 1 1 -1^X000 


ft 


CC HU.M -^,0 *— . » 


ft x xxxxxxx 


1 O 1 








LL 


1 


< ZtXftZXIXU 0 


ft 


— * *00— o^» 


ft h-<XCQ< JLL ^ C3 


1 O 1 


M 


O 






1 


0 1 ^00 ZZZU.QUHM2.X 


ft 


»— <1^^0U^00 0 


ft 


c_?or 00 


1 -J 1 




M 




LU 


1 


LU 2 r<T 2 . 


ft 


•-» OOLO — OLOO 


ft 


3 


1 O 1 




LU 




X 


1 


Z>>>> hhmmm>»>>> 


ft 


OvfOLn— oiriw * 


ft 


DC 


1 X 1 




CC 




H 


1 


O0O0O0O0O 00 0O0O0O0O 00 00 00 00 00 


li- 


*-ift *— <x— 00 


ft 




1 1 










1 




ft 


JJOOOQHZJ 


ft 




1 X 1 










1 




li- 


Q. *— 1 UJ <t O 31 


ft 




1 X 1 










1 




ft 


xlu— ajr>a>ic 


ft 




1 1— 1 










1 




ft 


o:ll dc 0 00 -joo 0 


ft 




1 1 










1 




ft 




ft 




1 I 1 




ft ft ft ft ft * ’ftftftftftftftftftftftftftft 


ft ft 




ftftftftftftftftftft-lt-ftft 


ft ft 



ooooo oooooooooooooooooo 



oooooo ooooo ooooo 



117 



21 - INVERSE Of Z 

2T1- INVERSE OF TRANPOSE OF Z 

X - AUX TRANSFORMATION MATRIX FOR COMPLEX EIGENVALUE 

XI - INVERSE OF LX 



ooooooooooooooooooooooooo 
oo o ^ c\j ro ^ m so co a> o r-4 cnj ^ vj- m no oo o o cm m 

r^ujujajcoajujcuuuujcuU'U'U'U'U , ij'U'U'<j^u>L^L_JCjc-J 

auyuoauQUuuugouoooouoHHHH 

ouoooooooooooooqooooooooo 

Z2ZZZZZZZZ2ZZZZZ2ZZZZZZZZ 

oaaauoauuuQOQaouoooouauucj 

Ui ^ >• ™ 1U( ILJl U lr M w n 



I 1 It 



* 

rM 



LL 

a 

UJ 

00 

cC 

UJ 

> 



uu 



N 



(03 

XLL 

q: h-ui 

CO 00> 

o 

(O 



* 



3 # CO 
CL-* 

5 !<h 
O I 3 
O’— 1 CO 
00 — 

X — lu<n 
►"Me > v s 

ocoo~* 

h- II CO | 

<X — < — 

5100 * 

— 00 * 
xococ 
uj * 

— JU I LU— ** 
clq:Zut 
Xlu <r — 

OX coco 

O j£ 

i i i i i i i a 

I UJ 

oocoocaca: on 



OOh- 

zu 

OUJ 
03 
— X3 

or cmoo-o 

OOQC X( 0 < 
a'O '-‘OcrTUJ 
coo I OoCX 

CL -fr — QC<h- 

co^o* + 

zz 

joo:m<m 



00 


iuioaotxij 


1 LJ 


CL 


(OCX 


OO -' JOQC 


t LL 


O 


00 


O co a. ui uj 


1 




CL 


OC ►— »H— CO 


1 UJ 




CO 


XUJ 


1 X 



QO 



O 



» 

oooo 

XX 

hh 

i— »i— i 

-S^ 







1 cO | 






QC 


1 ' 1 


OQ 




UJ 


1 00 1 


LLIUJ 




rM 


1 > 1 


H-U 




1— 11— t 


t 3 1 


^ oca: 




OX 


1 < 1 


— »>-ioa aroo 




a:— 


1 X 1 


UO<<QCOl“h 




*— *i— 


1 < 1 


DQ»OhOO 


00 


XQ. 


1 1 


00 00 00 00 h- CJ U 1 UJ — » 


— 


a 


1 UJ l 


00 00 — W(JUI»U 


co 


LL 


1 3 1 


<< UJ> — 




Oui 


1 3 | 


o ca:> ocuc 


•— i 


X 


1 <T 1 


aca^oo QC<<X 


h- 


h-l- 


1 > 1 


oohho:<jjo 


X 


UJ 


1 1 


hh0U<30D*H 


<CQ 


QU 


I a: | X 


UOUJUUOUUh 


— ) 


a 


i < i a 


UJUI»3UZZU 


CL # 


UJ 


1 3 I ►-» 


>> UZmmZ 




XX 


13 1 h~ 


QC0CX^h0000 3 




1— LJ 


I CO i o< 


a. a: <z <£ m oo ll 


CO 


i—i 


| Z UU CO 


<T <C -j 100 LULU 


1 c/i 


LLh- 


I »-H I UO — > 


3333 UI>>lL 


UJw I 


co 


1 00 1 # 3 3UJ 


3 3 o CO UJ > — M—» LU 


1— # — 


X 


| | CD 21 QC 


COCOXX>»— h-l— U 



COLL 
* # 
U-LL 



X 

a 



3 

UJ LLO ^Zmm,-h|-<<U 0 LL 
I JlOh H-.^oOOOh-<UOZ 
<QO c J Xoooo <(j^ M < aZ 

CL U J LLJ X *— • UJ UJ O ' O—JQC UJ 

ul uj 30 >— • qcujuj:>>-«— jclclh- ll 

OQ 3 JZ(Zh»MMjaMH OO 

UJ— KUJh- <»-"-• h- h- CL *-»h-K- h- X CO LL : 

UJa:<>0<Xh-l— I— » I — *»— 1 1 — — I— IXCOC LL# # zzz 
3 *— ' >■ m 2 ; mmQOH J 33 <^- c£ # UL CO •— m— ' 

—100 OUJ >CQQ 0 J 025 :JLLh CO — — 
<LliQU.I >hOO<OZ CL > > + + 

>ouj(Yghm<< s: t-h- +»— + x x 

X QCUiUJ*— 0 — I — | — h 30 ^ XU, mh i 

LJJh-UJOa:i— XI— h 03 HDQ.a — .< # 

C 0 3 OcXU» 3 LU 3 3 CLCL 3 CLI— I— LL —ICO 
»-*Q-a:OOUJQCLCLl~|-CLX 33 I CL I | | | I 

LUXO QCDM-Z 33 ZHOC I I 

I— • | O'—* mmCQm cl I 

It | 111 X XLLXH 3 LL 

I Oil I I I I I JllZLQ-Ljaj 

O 3 *— • »-• MMOOa,CL<hQ.h-CLiJia 

hOBUN^XhhUQS^Z^XhJJhJXhm 

tu 5 iaaox^o> 3 >o> 3 >a.xaaxcuLxx 






OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 



118 



ooooooooo o 


O 




OOOO 




O 


O 


O 


OOOOOOOO 


^mor s -oochO'HOsi r*. 


>r 




Nfr^coo 








<■ 


cm co m or-* coo' 


UUUUUUhhh «-h 










*— i 


O 


— 


ouuuuggu 


HHrH 


»hhhhhh h 


»*H 




pHsUnOnO 




-1 


O 


H 


OOOOOOOO 


ooooooooo O 


O 




0000 




O 


O 


O 


OOOOOOOO 


2ZZZZZZZ2 X 






ZhhH 




X 


X 


~Z. 


^zz 


zzzzz 


UOUUOQQQU U 


O 




O CL CL CL 




0 


0 


O 


uoouoaaa 


UUUUUUUUU L J 


O 




OOOO 




0 


0 


O 


OOOOOOOO 




OhhX 




















41 * 


— 1 — •'LL 


4f 


4V 4f 4V 4V 


41 41 41 


41 41 


1 41 


41 






41 


— CLOCL - 


4h 




4V 


1 


41 


1 


41 






tt 


XX^Q--* 


4f 


1 


4V 


Cl 1 


41 


1 


41 






41 


<LL— -»0 


4V 


1 


4V 


00 CL | 


41 


1 


41 






-ft- 


-aox^ 


4f 


1 


4V 


<h~ 1 


-*-*41 


1 


41 






-ft- 


— «.»— »;> * * — * 


4h 


1 


4V 


CL | 


CM (M 41 


1 


41 






41 


ox *— *o o 


4h 


1 


4V 


>0 | 


- -41 


1 


41 






41 


fH 


4h 


1 


4V 


-1 1 


CM CM 41 


1 


41 






41 


••^Oh w *0 


4h 


1 


4V - 


H-5: 1 


— 41 


1 


41 






41 


OOH *— —*0 


4h 


1 


4V — 


00 1 


— "-h* 


1 


41 






41 


*0 UL CD r— < r— ^ •» 


4h 


1 


4VO 


<a: 1 


QCqC 41 


1 


41 






41 


— * *0 H-Jh -— o 


4h 


1 


41 *— • 


Xll 1 


21X41 


1 


41 






41 


NO^-^d ►OhH 


4* 


1 


4V - 


UJ | 


00 41 


1 


41 






41 


OrH'-X^O'HNw 


4c 


1 


4VO 


X | 


QO 41 


1 


41 






41 


*— a ul x *-* — a »-• 


4V 


1 


4V -h 


Q-J I 


- -41 


1 


41 






41 


— X *-J 


4V 


1 


4V — 


U !J— 1 


——41 


1 


41 






41 


OX -J^XQ-ZN 


4f 


1 


4V-4 


a:o l 


CM CM 41 


1 


41 






41 


m ^dooz 3 - 


4f 


1 


41- QQ 


<ui 1 


- *41 


1 


41 






41 


— -*Q . 


4h 


1 


4V •* 


Jd 1 


<MM4> 


1 


41 




X 


41 


OQH-* “ZHUJO 


4f 


1 


4V — 


00 c | 


— 41 


1 


41 




O 


* 


*— tf-H ►OOd *Q r-H 


4V 


1 


4VO 


LUO | 


a: a: 41 


1 


41 




00 


*1 


1 1 1 —CD *—^ * — 


4V 


1 


4V-^ 


00 1 


21X41 


1 


41 




< 


41 


OOH *-^— *0 ■*** o 


4V 


1 


4^ - 


1 


0041 


1 


41 




M 


41 


*<-H ^-*0 *— 'O >-H o ^ 


4V 


1 


41-0 


LULU | 


00 41 


1 


41 


1— 


O 


41 


«J i— « *"H 


4f 


1 


4V-H 


CG—J | 


* *41 


1 


41 


X 


* 


41 


Oh— X — CL *0 *•— » 


4h 


1 


4V — — 


CD I 


— — 41 


1 


41 


*— • 


CM 3C 


41 


iruxaiiOHox 


4V 


1 


41 OO 


H-< 1 


CM sf 41 


1 


41 


OC 


X Q 


41 


CL » 5~ ►— « r— < r-H 1 


4f 


1 


41 XfNJ 


Ul*-i 1 


* •• 41 


1 


41 


CL 


CC * 


41 


3x-*>x— ac — * 


4f 


1 


41 - - 


Dd 1 


CM^41 


1 


41 


HH 


2 00 


■fr 


2.1LU - 


>V 


1 


•x —0 


1 


— — 41 


1 


*■ 


-X 


rvj 


4c 


o jooo 


4h 


1 


410fM 


> 1 


3041 


1 


41 


X 3 


CD 


41 


*— % *OOOOdH 


4V 


1 


41 


0 | 


XX41 


1 


41 


<00 


- O 


41 


* *o0 — 


4V 


1 


41 


IUO0 | 


3041 


1 


41 


2I< 


2 : a 


4i 


O— 4 - •»— %— ^0.0 


4V 


1 


41 OX 


000 | 


00 41 


1 


41 


CO 


3 


41 


in *— ooooo—* 


4V 


1 


41 — < — 


CD< 1 


* *4f 


1 


41 


MOO 


00 < 


•ft 


w Q r»< r— < r— 4 *W 


4V 


1 


41— — Nf 


<CL 1 


— 41 


1 


41 


•— » * 


< a 


41 


* *— %X 


4V 


1 


41 >0 - 


M | 


CM vt 41 


1 


41 


-X 


0 rs - 


41 


Z^JQihoOcOJ 


4f 


1 


41 <rv)CM 


ma 1 


- -41 


1 


41 


X3 


- O0-* j£ 


41 


• XdZhHHv - 


4V 


1 


41 — — 


<i- 1 


CMnT 41 


1 


41 


<00 


CM *— 4 M 


41 


—* U- 


4f 


1 


41-^00 


> 


— 41 


1 


41 


XoO 


2! O0—^ * 


41 


O *— * -clooqc cc o 


4V 


1 


41 OOsJO 


00 


0041 


1 


41 


<< 


cc oo>o< 


tf 


H^O^XCLOLUrH 


4h 


1 


41 -H — O 


><UJ 


XX 41 




41 


OoO 


XX *Q -M 


41 


-OHOdOOh - 


4V 


1 


41 -tU. 


2 : X 


3041 


j 


41 


* 


Q < r— i *<X * 


41 


0~< **-< * -O lUQf\j 


4V 




41 oa: - 


2.0— i 


0041 




41 


*<M 


3 *XO— * 0 — * 


41 


r-H -O — -^aCOr-lCL 


4V 


OO 


41 — < 


d^h 


- -41 


00 


41 


XX 


1—1 — h-«O0"H * —4 


41 


— 0*-«000 - — x 


4V 


LU 


41—0 


O O 


-*—*41 




41 


<a: 


LJCJ * 


41 


M^H*— X fHfH-^-^CNJU.1 


4V 


O 


41 ^O-H 


510 


CM CM 41 


LU 


41 


iZ3QQXX*-i> 


41 


' - -NOXh 


4V 


GO 


41 - 


OLUOl 


- -41 


-J 


41 


<OH-MQOOOQ — M 


41 


— X>‘-«00 — -HCL * 


4V 


< 


41 - -O 


ZUCD 


CM>t 41 


CQ 


41 


QoOO 


- X - > - 


41 


O CD — X 'HHh - •— X *— < 


4V 




41 — O-h 


— coo 


— 41 


< 


41 


^QUZh^ZO- 


41 


in *-* >— LUbClUCL 


4V 


Od 


41 O—^ — 


^000 


CQCQ41 




41 


*<MOhO *<— <Q * ^ 


41 


— *0 ►hXOQ / h- X 


4V 


< 


41 *— <— « 


Od 


XX 41 


on 


41 


XX * 


O ►— 1 4 


41 


OOH-.X3C -O *LU 


4V 


> 


41 -<*-* 


OCLOC 


3341 


< 


41 


<cdonooo>LO-<x 


41 


>-ir-4 *OLL — h“ 


4h 




41 OCCD 


-J LU 


0041 


> 


41 


ZZh-N^ZoO^-N 


41 


1 1 1 *OHh* fc O *0 


4V 


OO 


41 -H - - 


oaiz 


- -41 




41 


<000 


0 • *0 *^ 


41 


COO*— < — — J — * *"H *— * — A Z 


4V 


> 


41 — — » — « 


U.XZ 


— 41 




41 


00 *QnQLU — IdJXvO 


41 


41 -4wQdO -O -M 


4V 


l/> 


41 <roo 


h— •— i 


cm <r 41 


0 


41 


^HU--fCdZHri:UH 


•Hr 


H fcfHO HU 


4h 


r— 


*#1 aj >-H 


LU 


•* -41 


< 


-ft- 


0 


-ft- GL —CL -}h 


-ft- 


XX3 — -* --H *-hX 


4V 


CL 


41 - - 


xzu 


^N<r 41 


CL 


41 


QCUXX OCX XX 


41 


IUCO -MOO^OwLU 


4c 


O 


41 TOO 


h— •— «h- 


— 41 


X 


41 


LLJOOLUli JLU LUOOUJIUU1 


41 


_J *— * ^hhuThm J 


4V 




41 41 *-H-H 




<<41 


*- « 


41 


O O 


OOO OOO 


41 


CL— 03 *— O — *-CL 


4V 


1 


41 -J — — 


1 


XX 41 


0 


41 


luoclclcqluocgclcl 


41 


2.0*-i -OX^C/MZ 


4V 


| 


41 <OZ 


1 


33 41 


1 


41 


H-CXX3h-<3XX 


4* 


O —< *— * r-l ^ »-H O O O 


4V 


1 


41 UJttTO 


1 


OQ41 


1 


41 


ziuooozujnoo 


41 


O -00—0X3 -O 


4V 


1 


41 CTlCL - 


1 


- *41 


1 


41 


— alOOO — uCQOO 


41 


• h c\j m vf 1 n v-o r*- co 


4V 


1 


41 'Hcg 


41 


n n* 


1 


41 






-ft- -ft- 




4h 4V 4h 4h 4V 




41 41 


41 41 


-ft- 






OOOO 


OOOOOOO < 


OOOOOO OOOOOO 


1 



119 



o 


ooooooooo ooo oooooooooo 


o 


oooooo 


O 


oo 


OOOOOOOOO 


<r 


o co o ^ c\jm vt m so co a> o ^ og m vT to so r- co a> 


o 




CO 


OnO 


cm ro vj- m no r- co 




^^l\4iN^lN4CNJiN4lNiNi\i<V4 n|n|^uxir»i;M^in | 


i^int 


NT 


sj-vj- v^-vT 


>r 


sj-li ) 








'■H'-ft »“H <— < «— 4 <-H «-H H <-H *"M <“H •“H >— 


-H *H 


1 






H 




H 




H H H H H 


<— < P-M r-H «™H 


o 


ooo ooooooooo oooooooooo 


a 


oooooo 


o 


oo 


OOOOOOOOO 


X 




22ZZX2XZ2ZZZXX2Z2X2 


X 


xxxxxx 


X 


XX 


ZZXZXX2ZX 


a 


uou uouauauuoauouaaouua 


u 


oooooo 


o 


oa 


uaaaauQuu 


o 


ooo ooooooooo oooooooooo 


LJ 


OOOOOO 


o 


OO 


OJOJOJOOJOJUOJLJ 


ft ft O0ft ft 




o 




■ft -ft 


ft ft 


ft ft 


ft ftft ft 


ftftftftftftftftft 


si- 


x i 




Si- 


in 




-ft 




ft 






00 H- 




ft 


a i 




ft 






■ft 




ft 






h- a 




ft 


►—» 1 


< 


Si- 


< 




-ft 




ft 






a -j 




si- 


h— 1 


h- 


ft 


X 




•ft 




ft 






-j a. 




ft 


o | 


< 


Sfr 


•» 




-ft 




ft 






CL LL 


OO 


tt 


3 1 


a 


Si- 


— * 




-ft 


1 


ft 






UJ 


*— 1 


si- 


3 1 




ft 


o 




-ft 


1 


ft 






UJ 30 




ft 


LL 1 


H 


Sfr 


in 




-ft 


1 


ft 




• r 


3 JX 


Q- 


* 


1 


O 


■ft- 






-ft 


1 


ft 




00 


—1 <.►-• 


UJ 


ft 


CL | 


CL 


-ft- 


X — 




•ft 


1 


ft 




< 


< >00 


h— 


tt 


Ul I 


X 


■»■ 


QO 




■ft 


1 


ft 






> 3 


00 


ft 


X 1 


»— 4 


-ft- 






■ft 


1 


ft 




00 


a: 




tt 


h- 1 










tt 


1 


ft 




LU 


aC <lu 


OO 


Si- 


O 1 


X 


-ft- 


Oct: 




tt 


1 


ft 




Q 


<LUQ3 


*— 4 


ft 


1 


o 


■ft* 


tn - 




tt 


X 


ft 




O 


JOZJ 


X 


Si- 


O I 


CL 








tt 


LU 


ft 




X 


3-JcO< 


h— 


ft 


X 1 


u. 


■ft- 


oo 




tt 


H 


— » ft 






o<x> 




ft 


< 1 




■ft 






tt 


00 


^ft 




> 


X>*-^ 


•> 


ft 


1 


00 


■ft 






tt 


> 


— ft 




LL 


co cL 




ft 


00 | 


□C UJ 


# 


— > 




tt 


00 


<ft 




►— i 


S)cL < 


LU 


ft 


Cl | 


LU -J 


■ft 


O LU 


X 


tt 




— ft 




h— 


<X-J 


a 


Si- 


LU | 


fxJ CO 


-ft 




o 


tt 


UJ 


Qft 




X 


O-Jl-3 


< 


ft 


H 1 


-« < 


•ft 


*—« •* 


o 


tt 


X 


II ft 




UJ 


>3^0 


Q O 


ft 


LU | 


X — 


■ft 


o— » 


o 


tt 


h— 


Oft 




a 


^O^X 


IU UJ 


St 


X 1 


l-H IX. 




QO 


o 


tt 




ft 




►— • 


X ~ 


oo a 


Sfr 


< 1 


h- < 


■ft 


•—•in 


\ 


tt 


CL 


•ft 






00 ►-i 00 CO 


3 


tt 


CL 1 


CL UJ > 








tt 


3 


UJft 




a 


OOOOOO 


1 — 


* 


< 1 


OQQ a. 


tt 


—a: 




tt 




•ft 




h- 


< <LU 


UJ 00 


*«■ 


CL | 


XX >- X 


tt 


CJ— • 


*— « 


■ft 


h- 


*-«ft 






cl UQ-h 


CQ CL 


Si- 


1 


LUDGQCh- CL O 


tt 


in - 


QC 


ft 


LU 


ft 




— Q 


X 3 


—4 •— • 


ft 


<X 1 


IOOO o ~ 


tt 




X 


ft 


00 


ft 




muj 


(-31-a 


OO UL 


Si- 


OJ | 


HaSCQhH I— > oo 


tt 


03 




ft 




CLtt 




•-4O0 


00 cOX 


• 1— 


Si- 


r -sl 1 


OX 0< LU 


■ft 


► H 


> 


ft 


a 


ft 




•*3 


CLLLJCLC 


O CL 


ft 


M | 


a^L^ujM lu ql a 


-ft 


«k 


o 


ft 


H- 


+ tt 




in 




H> O 


Sc 


X 1 


QlULUX > jCUJLU UJ 


tt 


OO 


— o 


ft 




Xft 




*—*00 


LL — UL 


UL LJ > LL 


Si 


>— 1 1 


LL "2 CL QCOC <OOOJUlX 


tt 


m-H 


oo 


ft 




LJLft 




— " OC 


X 


mxo 


Sfr 


1- 1 


O CL l — fr— * LU H— <<00X*~*UJ 


tt 


• — ■ •— * 


— »o 


ft 


< 


1 ft 


o 


UJ 


X^— * Xh— 


-'LUX Q- 


ft 


Q_ | 


UJ _J 3 X oo 1— x LU CL CL O •-« h- X 


tt 


CQLL 




ft 


CL 


II ft 


*— * 


H- CQ 


0H02 


hJLULU 


ft 


a i 


O ^XLULUOOOOOOH-X-^ 


tt 


3*-* 


h-X 


ft 


O 


3ft 


\— 


<X 


a. lu 


<OJh 


Sc 


i 


CQU1UJ<aX*-< X30l— 


tt 


> - 


> X 


ft 


a 


ft 


< 


X3 


CLQOCX 


XUJCroo 


Sr 


QC | 


<OOiXO<QnCX XOQOC3 


tt 




cl ;i 


ft 


CL 


QTft 


oc 


OCX 


Ul LULU 


QCQ^LU 


ft 


O | 


•— «lgcq i — oc<OQCa: X O 


tt 


— O 


^-*ac 


ft 


Q_ 


Oft 


< 


o 


*-*x*-*o 


OLLQC> 


ft 


LL | 


0C<X<I00X<C o^hQOU 1— 0C 


tt 


OO 


•^X 


ft 


X 


ft 


CL 


LL H— 


XLUX< 


LL LLO 


ft 


1 




tt 


mo 


— • » ^ 


ft 


LUO 


ft 


UJ 


3 


— x-^-j 


X X 


Si- 


OO | 


->cf. o:q*-« fCUH 


tt 




ft 


X 1! 


ft 


QC 


•*CL 


1 — 1 — h— Q_ 


• -3c JLU 


Sfr 


-Z | 


<<oi-ixa:ccQCQCOuooai 


tt 


CD — 


HO 


ft 


h— > QLft 


CL 


J2 


CL CL 


-*XX3 


Si 


a i 


X>> <TLU UJU )UJ U : tt 1 00 o 


•ft 


-JX 


•»tn 


ft 




ftft 




ct:*-* 


OtUOUJ 


— 1 •— •*— *Of 


Sfr 


t— 5 j 


O cervix Mrvjrv) >00 


•ft 


>^os 


ft 


O •• 


‘>ft 


h— 


h* 


h- -J 


LU X H— UJ 


ft 


00 1 


^ZZUI^ihmm^mDVCQ 


-ft 


►*— i » 


ft 


h- 


OLft 


< 


XLU 


a:3a:a 


Q <T QL CL 


Sfr 


X 1 


LOOONrtrrss cq > 


■ft 


«k 


—o 


ft 


X3 


1 ft 


X 


Ocx: 


UJ CLUJ CL 


^x<u. 


ft 


LU | 


LU*— •►— * *— i »— 1»— i *-h >—<*—• h-h>. QcQ 


-ft 


o— 


zz 


ft 


>-<CQ II •»• 


ac 


x< 


zozo 


♦-4 h- 


ft 


X | 


QoOcO X I — h- h~ 1 — 1 — h-^HI LJ 


■ft 


lOqX •- 


ft 


+ 341- 


o 




21^1 


Si- 


►— i | 


LULU •— 4 Q» CL CL Q clclqllu^jlj 


-ft 


o*-^xx-ft 


<IX 


ft 


U. 


-OOLUOUJX 


•» 


ft 


a i 


luqoh-ooooooxooxlu 


■ft 


xo 


•^OO-ft 


l-<L Lft 




UJqCOC 


XX l 3: 


ft 


i 


X CL <3 00 


ft 


o 


— • m i-h -ft 


<X II QCft 


< 


Qhh 


< C 


Si- 


CL | 


1 — t 1 C3 1 1 1 t 1 1 13 


tt 


X— »0>sX-ft 


LJ^UJ* 


h— 


□ 22 


X X— *UJ 


ft 


O | 


1 1 


tt 


LJ— <* 


-H- 


O -L ft 


<X 


^UuO-hh^ 


JLS.L J 


Sfr 


1 


l com id -hX >- 1 


tt 


»— iH 


-cxo-ft 


1— 0_ift 


O 






Si- 


h* 1 


03 QOUL<^XLLc^LU 


tt 


LOXOXO^ 


3X ft 








rg 


ft 


LU | 


X»OhhQX^hhm^ CL 


tt 


X^r-. 


ft 


CL 


ft 










Si- 


00 1 




tt 


LU •> 




X 


ft 




ULULLI 


CJ 


ft 


1 




tt 




Sl—H-ft 


*— 4 


ft 




ocaa 


a: 


Si- 


1 




tt 


— *00<<ft 




■»<- 




<QOOO(\Jrsj 


< 


ft 


1 




tt 


am 


•— • Q Q -ft 




ft 




02zz: 


o 


Sfr 


1 




tt 


-H<\| 


tt 




ft 










ft ft 


-ft- ft -ft- 


ftftftftftftftftftftftftftftft 


tt 






ft 




ft 


ft 


ftftftftftftftftftft-IVft 














ft 




ft 











' ft ft 

■it 

UUUOUUUUUUUUOUUUUOUUOU OOO OOOOOO OOOOOOOOOO OOO 



120 



ooooooo 
o morgen in eg 
suunsu^J 'uuu^- 



oo ooooooo ooooooo ooooooo 

m in o co in o ua '-h rg m vf tn o co o o tn ^ c\i 
p» r>. p*. p*. p* p» u \p- cu u \au cu uj lu oo co uj lu au lt* u ^ u> 

^hhhhhhhhhhhhhhhhhhhhhh 



oo 

CO<f 
O' O' 



OOOOOOO 

ZZ22ZZ2 

uuuaooa 

ooooooo 



oo ooooooo oooooooooooooo 
uoooooooooooaa ooooooo oa 



oo 

zz 

00 

00 



0000000000000 
inm^i^coa>o~*cgmu^Nr-n 
u^i/u^u^uuDuinuu 
—i rH -H — 1 <-1 f\J < M CNJ 1 M AJ C\J 
OOOOOOOOOOOOO 
ZZZZZZZZZZZZZ 
uoouauQouuauo 

OOOOOOOOOOOOO 










H- 














Z 

UJ 








OC 














Z 








LL< 


O 


















X 


OX 


LU 












xa: 






a 


O UJ 


cC 






or 






OUJ< 






(—4 


• X 


UJ 




CO 




LU 




h-<tCL|- 






tOH- 


>UJ h- 


Q 




O 


O 


X 




COLU C 






LU<t 


• co 


►— » 




<r v) 


CC 


a •— 




3 toa 






-JO 


C0< CC 


CO 




a 


z 


Z K 




200 




h- 


OO 


X 0 


Z 




UJ < 




< 




h-aCLU 




LL 


CLX 


UJCLX CCCC 


O 




UJ 


LL 


< 




X <x 




»— * 


a: 


X —* CC UJ 


O 


0 


CO UJ 


-J 


co 




*-HQfOh* 




X 


LUO 


l—OO LUN 




LL 


LU 


O 


LU H- 




cLlu 




00 


OIL 


ZaC 


h— 


Z 


a. 0 co co 


O 


O < 




h-u J 






Cl 


cC<< h-X 


Z 


H- 1 


X < 0 


z 


*— 1 h— 




< <LUCX 




UJ 


•-*»- 


— ox Z— * 


M 


•» 


O 3 LU < 


•* 


QfZLU 




Xaroci- 




a 


h-Z 


in uuz ujh- 


< 


1- 


O X 3 


LL 


h-LUZ 




LL HJm 




< 


— J LU 


H-I XQ- 


CC 


z 


< a uj 


3 


<20 




QHo>^ 




0 


3t0 


•— 0<i-4 LU Q 


H 


(—1 


l~ X ^ UJ 


a 


Xlu 




XOLU 




UJ 


51 UJ 


O UJO< O 


co 


a: 


Z Z CO 


OC 


OLU 




ZLU»-4C0>- 




a 


□C 


• X> O <LU 


Z 


CL 


UJ CO X 


z 


1 — <t O 




azz 0 






LU CL 


OCLUO OX 


O 


H-t 


*-h 0 coo co a: 


* 


300 




a: lu i— 




X 


X 


HQ^ J<QmQ.H 


O 


•* 


0 <OLUOUJ 


0 


CL CL CL 




LUh- CQO 




0 


H-LU 


ILO COZZ 




0 


<C0 <th— <h" 


0 


z 




XXZ UJ 




< 


X 


coOLUa:<^-Luo 


UJ 


UJ 


CC IU LU IJJ 


u 


hhcLLU 




oi—Lu^a: 




UJ 


Zh 


— a:ouj Xoh- 


CL 


z 


Oo LU LU X U 1 X 


0 


OO 




< x<»— 






f-i 


b-a--J:>uj>oz 


>• UJ 


0 


CO COUJOLU< 


z 


LLLL < 




UJ <XXQ 




h- 


z 


<1 <1 •— ‘COCOCL*— • 


h— CO 




0<ILO co QTCO UC 


•» 


u 0 




O 




<2 


Q ►— » 


XcO>Z< 


O 


*1 


Z»-*h-> <2 <C 


0 


|— CL 




QCQCOUJO 






— LU 


CHhh DlUKh 


CO 


h" 


M ^ZOQtU Jd 


z 


COO 


• • 


OO a:z 




x 


OtOQ 


ax<x clozo 


h-a: 


CL 


Z<*-'LLJOc0O 


0 


Z CL O 


CO 


ULLLCOLUO 




< 


• OuU 


LL 1 — <IO UUCL 


LUO 


0 


QC Xhl CC O 


a: 


rzh 


< 


—•XCL 


z z 


a: 


O cO 


oxozo 


COLL 




LU LL*< 1— Oh— O 


z 


3 1 — • 




cOQX) — co 


00 


O 


*-HLU 0 


••ZD^LUOaro 




•» 


OZhdiiJ<Zd 


•» 


O CO 


CO 


Q<h LU 


a: a: 


O 


LL CD 


O'— 'OtOac;*-«UJh- 


h— O 


1— 


ZO(/)hZZOh 


CQ 


□ uiZ 


z 


~^^nC\u -qC 




a: 


rO h- 


0 ZO^cocl 


<< 


< 


OmZco Ooz 


O •» 


OLD < 


OCMCXIC OlLL Of Ql 


>> 


Q. 


-OO 


mQh COQC CO 


ID 


CL 


ucou Z^-H 0 


03 


LU 


CC 


- **o OLUO 


xao 




h-h-Z 


ZLUCOQUJUJHLU 


h-Z 


J— 


UJUCCOoOHO 


ox 


0 ax 


^c\ J UJ CQO 




> 


<r 


*</) 7A>UIX 


<r 


CO 


(XO uimZZ 


zz 


Zh 


> 


iuqCqL X • 


zz 


CQ 


XcO — 


ODQD coco 


OCX 


» 


LU LL M| — LUt — 4 QC 


•» •» 


< 0 


CO<<af<LULU|— UJ 


00 




oca: 


CC U_l O <1 1— *— 1 


LU 


•1 


HU.Om<2:c£LU 


CQcL 


CO 




O 02. O 


CCcC 


0 


OO cO 


•oOaTLL DQ5 


1 — cO 


z~ 


•UJO ZM-CLN 


zo 


touj II 


Q 


XCOZLUZ— ' 




LU 


LL h— UJ 


O'—"-' • U J CL H- 


LUO 


amx cc^mq h-. 


00 


Z-J 


<t 


Q— QfOLL 


OQ 


HH 


00 


Z LO — >h- CC 


x< 


O ' <1 QC LU H— X COX 


CC 0 


003 


LU 


a: cc<i o 


<< 


LL 


••<0 


«oujomdo:uj 


< 


ZODlX:UJCQQ_i— «LUO i~» 


zz 


OCCLX 


OC 


— u<x • h- 


LL'LU 




rOU-X 


X OO LL OQOZ 


OTLU 




• <CDXOl— Z<h- 




^ 2! 




OO OUJOcoO 


of a: 


O 


1- 


+ J— O 


< UJ 


>h- 


•CLXO CL O CL 


<d^U. LL 


• 


- hu cl of qccl 




O 


ZOh 


COZ | *-«0 1 O— * 


Q_CO 


o< 


oz 1 0 10 


O z »— *o O •• 




Hi\JZd <IUZ 




5! 


"ZZ 


•» — <t cO 




zx 


: 1 z 1 


OOM uj 


X 


— — OLuiuocn^ 


XX 




CNJh-iUJ 


HH >-NH CD LL LU 


1 


•*QC 1 H 1 H 1 


oof — 'XoC O 


HH 1 


<< clX X 


mm 




h-h-x 


ZJZ Z Q 




OOO 1 < QZ 


ZZ 


— LULU a. 


Of 


LUQh-<3lU 


cCcC 




-^XLU 


H^H 1— 1 1 


O 


<TLL<t za:hiuMa 


- ajcnx 


1 — 


coa: zx 






*00 


xox X 


O 


CC 


a:>OhdZ(vii 


<o<xx< 


< 


CL *C0 h- 


<< 




— ' »<I 


> > >0 


•— » 


O 


ooocoouaz 


ZOX30X 




O Z<ui 


XX 




h-ui-l 


to — 00 coa" 


z 


•— 4 


•— <ZZ»— ^ <•— 1 


oa^zziu 




OUJ*-» XZ 






ija. 










LL O U 


<L 


MX 1 | 


UJO 








1 






Z X LL 




£h-OUJ 






CO 




1 


in 








r- 


<coao 


00 Q> 








1 






nO 






LUCCLU co IJJ 






a 


0 


1 


0 








0 


CC LUDLC h- 


OQ 




a: 


cC 


1 


CL 




O 




a: 


LUI— h-H- CO 


CCCC 




< 


< 


1 


< 




cC 




< 


ILLLU<UJm 


<< 




0 


O 


1 


0 




< 




0 


i— < caoco 0 


00 


* -ft- 




* H 


1* ■H ■tt’ 






* * JV ^ ^ ^ ^ ■* 



121 



o 

NO 

LJ 

eg 

O 

2 

o 

o 



ooooooooo 

in 

C _J i—4 —4 f— 4 «— 4 f— 4 *— 4 

rgrgigt\jrg<g<N<NJcg 

OOOOOOOOO 



UUQUQUOUU 

uuuuuuouu 



ooooooooo 

^Of^^G0(>O^rv| 
ig —4 *— m g f-H >4 V g t. g 

<N eg i g ( m c\j c \j i g i g rg 

ooooooooo 

ZZZZZZZZ2 

□oouoaooa 

ooooooooo 



oooooo 

m>r f\J in Osl O 

egigtgcgigig 
eg c\J eg rg rg ( g 
OOOOOO 



OOOOOOOO 

r-coo'0^c\irrixr 

egigigiMfiMMiMin 

rgrgrgrgr\jrg;grg 

OOOOOOOO 



ZZZZZZ ZZZZZZZZ 

oooooo oooouuoo 

OOOOOO OOOOOOOO 



##**»#♦###**♦*♦##**#*#*##**#* 




Jf Jf 


Jf 




1 






UJ 






•H* 




Jf 


Jf 




1 


~*o0 




hi 










Jf 


Jf 


1 


1 


LULU 


UJ 


oof— — * 










Jf 


Jf 


1 


00 


20 


O 


CL X *H 










Jf 


Jf 


1 


f— 


— «0 


CO 


►—•OOO *»— « * 










Jf 


Jf 


1 


2 


h-< 


< 


u — j — • rgrg 






* 




Jf 


Jf 


1 


4—1 


LU “D> 


M M 


_JOO — r. 






* 




Jf 


Jf 


00 


< 


DhO 


2X 


UJUJUJ 0—4 










Jf 


Jf 


H 


CL 


o2clcd 


<2 


IhU LL ^ 




2 


*«• 




Jf 


Jf 


X 


H 


< LLI 2! 


>0 


H* U.CLO 




O 


•If 




Jf 


Jf 


CD 


00 


>Xcl— • 


LL 


o< 00 


00 


►—4 


# 




Jf 


Jf 


*—4 


2 


UJL’Jj— 


O 


— CL 00 


UJ 


H 


•If 




Jf 


Jf 


LU 


O 


2022 


O LU 


— 4UJ|— LU— LU 


O 


< 


•If 




Jf 


Jf 


2 


O 


<<2< 


-•X 


►MO 2rg2 


O 


0 


•If 




Jf 


Jf 






JjMh 


00 H- 


1—4 0 0 ► < 


CL 


0 


•If 




Jf 


Jf 


CL 


LL 


13 a oo 


UJ 


'-'uj —4 x: 




0 


•If 




If 


Jf 


LLJ 


O 


o 2 : 


OLL 


aiuo 


a 




■If 




Jf 


Jf 


NJ 




2:000 


O 


h-— • JCJ> 


Ul 


UJ 


4f 




Jf 


Jf 


•— 4 


LU 


MJ(XZ 


UJ 


UJ 2 0 


CL 


0 


•If 




If 


Jf 


X 


CL 


00OLU <— 


XLU 


2 •— o>o 


4— 4 


a 


•If 




If 


Jf 


t— 4 


> 


CL'-' 2: 


1—0 


< XLU 0< 


00 


CL 


•If 




If 


Jf 


K 


h- 


a: luo 


X) — 00 • *»2 O 


LU 




* 




If 


Jfr 


CL 




OUJU.XUJ 


OOQ 


CL< Q O — 


a 


a 


•If 




If 


Jf 


O 


O 


LLX jCOO 


UJO 


oauiuo h- 




UJ 


•If 




Jf 


Jf 




2 


f-2 O 


ZD 


1— LL CD— « •< 


X 


— CC 


■If 




Jf 


Jf 


•k 


< 


> uiO 


-.a 


00 IT CLLU 2 : 


0 


0—1 


-If 




If 


Jf 


a 




-JoO X II 00 


' U.3 


ccaH4<h «o 


< 


• 00 


■If 




Jf 


Jf 


UJ 




ZliJhQ^ 


LU 


ujooriuj^h 


LU 


OLL! 


■If 




If 


Jf 


cc 


O 


OO 


o>* 


M t— ox O 




-40 


Jf 




Jf 


Jf 


LL 


LU 


OO— 


CL 


^o x • < 


2 


LL 


•If 




Jf 


Jf 




00 


oxo 


Xh 


ZZZUJhQ 


O 


<N> 


■If 




Jf 


Jf 


< 


O 


3 LU <3 


3 ~2 


4— 4 <r *— ‘C0ODUJ5: 


LL 


— CL 


^f 




Jf 


Jf 


I— 




a^)ai a 


O CLUJ 


1— ^ 2< 


— 


<t 


■If 




Jf 


Jf 


O 


h- 


2 0UJ02 


23h- 


CL— OUJ »— »QC 


OO 


h-2 


■If 




Jf 


Jf 


UJ 


a 


X 


a<uu 


a — — H-aDOO 


• 2 


<H-4 


■If 




Jf 


Jf 


Q 


2 


>00h-LU> 


>2XX 


*— < —4 OO 


o< 


IO 


4f 




Jf 


Jf 


*» 




rQt-iOK-CQ 


CO h- 


LU— ' ►HOli.'cC 


H L) 


CL< 


■If 


Ch 


Jf 


Jf O 


> 




< 


>00 


XXNQhOCCL * 


LL 


0X0 


■If 


O 


Jf 


Jf QC 


O 


00 


3UJ223 


3 CQ 4— 4 • 


h-w^z: 


oguj 


LL— iQL 


If 


nj 


Jf 


Jf 1 — 


23 


LU 


OiXOUJO 


O la: 


x Oooiuo 


--2 


< 


4f 




Jf 


Jf 2 


LUO 


ID 


CL LULL 2: QC 


2-£h- UJ 


O LU 00UJ*-4Xh- 


1 — O 


•*QU 


•If 


O r— 4 <— 4 


Jf l/)Jf OO 


ZDium 


O 


X LU 


O M 


001 — —•cL 1 — UJ 


< 


2 


-If 


► 


Jf UJ Jf CL iL 


OQh 


< 


O *00 


02 ■ 


LUCQ<I ►—X 


2TUJ 


=)<a: 


Jf 


O 


Jf O Jf h* •» 


UJ 


> 


<QUI < 


c X 


H* < X CD O 2 • — 


OCX 


X UJ 


4f 


1 1 1 


Jf Olf 2- 


cLO> •* 




LLUJGJUJ 


LU • *h-4 


MhzaiuQH 


0 


CD — 1 CL 


■If 




Jf X Jf 00 


LL— »rg 


QC 


22 002 


2 CNJH- 


QcL <UI|— LU ► 


LLO 


<< 


■If 


sQ gj 0 


Jf 


Jf 20 


3h- 


< 


LUZ^h 


X — CL 


o<cL 2 20 g 


O 


X U IUJ 


■If 


in in m 


Jf h- Jf ►CL 


O *3: 


O 


3** 


• ► 1—4 co CD 


O > UJ LU OO—' 


• »»--4 


CL— 1 


Jf 


csjrgnj 


If ID If UJ h— 


<x - 


O 


xzuo X 


X 2 ■*-' 


0 N J00<mO 


cq:s 


- a 


■If 




Jf a if 02 


•— 4<t r-l 


CD 


•— ‘ILJ ►— 1LU*— • 


*-4»— KUJfMcr 


> • zoo 


X 


IDl — D- 


■If 




If 2 If OO 


h-Xf- 


2 


CL 002CL 


2<<X 


oz< a ooo 


>u« 


XO 


■If 


O CO 


Jf — 


* Jf 2 2 


•— 43 : 3 


♦— 4 


KUL>-— •!- 


hxri-HO^MMMao« 


**CL 


OCLLU 


-If 


PJCNJ— 4 


Jf 


Jf - 


2 


00 


< _jH<r 


< CL 


— iO0h-2 < 0 


XU l 


<22 


Jf 


'-W 


Jf 


Jf — 






ZUJ<DZ 


xrsoo 


XLU CL < OOOOOO 


ox 


UJH-fO 


Jf 


1 — 1 — h“ 


Jf 


Jf O — 


OO 


0 


IZO 


02*- 


HOOXOmUJ 




CL 


If 


UJU1LU 


Jf 


if injf 


•— ih- 


LU 


LL 1 — <20 


2— • 








If 


ootyju) 


Jf 


jf p- ► 


ZMI^ 


CL 












If 


LL QLUL 


Jf 


►<? 


X •* •" 


t—i 


O *-4 


<\im 






in 


Jf 


qc cl cl 


Jf 


jf — 4 —' 


<r^4 H 


OO 


—4 —4 


>— “4 r— f 




rH 


1—4 


If 


LLi U *LU 


Jf 


Jf — 


X — — 


UJ 












Jf 




Jf 


Jf LU 




0 


O O 


OO 




O 


a 


If 


OOO 


Ji- 


Jf Qh 


f-oa 




2 2 


22 




CC 


CL 


Jf 


ooo 


Jf 


Jf <— • 


o<<^ 


h- 


< < 


<<r 




< 


< 


Jf 


<<< 


Jf 


Jf UJCl' 


CLLLUU 


X) 


O O 


00 




0 


O 


Jf 


ooo 


Jf 


Jf CL 2 


22 CL 


CL 












Jf 




If 


If 


►—4 


2 


Jf Jf Jf Jf Jf Jf 


Jf jf jf jf jf Jf 


JfJfJfJfJfJfJfJfJfJfJfJfJfJfJfJfJf 




Jf Jf Jf 


1 


4—4 



I I 

I I 



oooooo OOOOOOOOO oooooo ooooooooo 



ooooo ooo ooo 



122 



ooooooooooo 
in no h qo on o <m cn >t in 
>r 

cNCNCNiNcNjcgcgcgiNCNcg 

ooooooooooo 



□uouuuuoauu 

ooooooooooo 



oooooo 


ooo 


OO 


ooooo 


ooooo 


oooooo 


'O cxd vjN v0 r- * 


0~*<N 


m>r 


nsor-aoa> 




sor-aoo^<N 




uunut 


u \u \ 




r-p-r-p» 


rr-r-cuujuj 


c\jc\|c\ 4 c\jc\jc\g 


rgogrg 


V 4 CNI 


rgrgrgrgrg 


c\jrg\|04f\i 


CMOsjrgcMAjrg 


oooooo 


OOO 


oo 


ooooo 


ooooo 


oooooo 


2 Z 2 ZZZ 


zzz 


zz 


zzzzz 


zzzzz 


zzzzzz 


uaauaa 


ooo 


oo 


ooooo 


ooooo 


OOOOOO 


oooooo 


ooo 


oo 


ooooo 


ooooo 


oooooo 





i 

l 

i 

i 

i 




1 

1 

1 

1 


* 

LL 






1 

1 

1 

1 

1 




1 

1 

1 






1 

1 

1 

1 






i 




1 


-J 






1 




1 






1 






I 




1 


O 






1 




1 






1 






i 


o 


1 


o 






1 




1 






1 






I 


LL 


1 


X 






1 




1 






CL 






i 


X 


3 


•k 






1 




1 






H 






i 


HH 


cl 


LL 


1 




1 




1 






CL 






i 


•k 




3 


1 




1 




1 






00 






i 


1“ 


3 


a 


1 




1 




1 






> 






i 


Z 


or 


oC 


1 




1 




1 












i 


4—4 


•k 


X 


1 




1 




1 






LU 






i 


cc 


LL 


•k 


1 X 




1 




1 






X 






i 


CL 


•> 


a 


1 a 




1 




1 






1— 1 






i 


4—1 


O 


-j 


1 M 




1 




1 






H 






I 


•k 


«k 


a 


Xh 




1 




1 






X 






i 


o 


CO 


o 


* — << 




1 




1 






o 






l 


LU 


•k 


X 


QCO 




1 




1 






ac 








z 


< 


•* 


H— • 




1 




1 














o 




o 


<LL 




1 




1 






-J 






< 


►— « 


00 


X 


Xhh 




1 




1 






00 






x 




LU 


a 


o 




1 




1 






X 






zr 


h 


O 


a: 


aa 




1 




1 






4—4 






< 


CL 


t— ^ 


X 


X 




1 




1 












2! 


O 


a: 


•k 


LU 




1 




1 






o 








*— t 


h- 


aj 


X * 




1 




1 






X 






oo 


•* 


< 


-J 


HO 




1 




1 






*— 1 




o 


o 


1 — 


x 


o 


LU 




1 




1 






OO 




o 


< 


<T 




o 


z~o 


— * 


00 




1 






X 




4—4 




QC 


LU 


X 


H- lOC O 


o 


Q 






■— 








z 


LU 


H" 


X 


•» 


X<-J 


-J 


X 




00 


►— 1 


— • 


oo 




•k 


LU 


oo 


H 


CO 


H>C 0 


o 


X 




LU 


W 




UJ 




x 


00 


»— « 




3 


K-« 


u 


a 




X 


X 


w 


X 




a: 


— 


•k 


Q 


o 


3: LUO 


X 


CO 




—I 


X 


X 


-J 




•» 




X 


z: 


OCX 


cox 


•k 






< 


o 


X 


< 




a 


< 


a 


< 


zz * 


00 *— ( i 


r-4 


UJ 


— * 


> 


< 


o 


> 




z 


H 


o 




••X Of cL 


tuo 3 


II 


-J 


X 


X 


X 


< 


X 






< 




</) 


< — <r COO LL _J -J 


on- a 




CO 


W 


UJ 


1— 1 


X 


LU 




x 


Q 


•* 


LU 


JC^JJJJQQ 


CO -J 


•k 


< 


CO 


o 


•* 


►—4 


o 




> 




l> 


X 


□jooaaoo 


<00 0 




t— « 


X 


t — 4 


— * 


•k 


4—4 




00 


h- 


Q 


-J 


ooooooxx 


»— <4— » a 




CL 


> 


LU 


►—4 


— • 


UJ 




•k 


X 


z 


< 


zuzzzz - * 


OC X 




<1 


•k 




W 


4—4 






4—4 


CL 




:> 


►2 »c g rv 


<0C 


X 


> 


— 


o 


X 


W 


H 




Z 


X 


Q 




< »<COOU. 3 ‘S 


> LU 


w 




X 


LU 


X 


X 


X 




4—4 


4— t 


< 


QC 


SOC^^jC^OO — 


QX 


3 


X 


W 


CL 




X 


CL 




2 ! 




a: 


LU 


osooaoQca: -> 


ZZh 


Of 


o 


CO 


t— < 


< 


-J 


X 




> 


cl 


o 


H 


^Ooc^QCarzzoca: - 


o< 


ar*~i 




-J 


oo 


LU 


< 


►—4 




00 


LU 


4— » 


LU 


zazzzz - - 3 _l— • 


h-» a: 


3 — 


00 


> 


UJ 


CL 


UJ 








M 




X 


x - - * -3300— 


00OO 


o 


UJ 


> 


a 


X 


XQC 


UJ 


X 


— * 


*—» 




<1 


^ajoLuaaoj; 


LU LL 


a:— 


o 


a— 




X — 


X w 


X 


X 


o 


X 


o 


QC 


ao— — -www zzq: 


ox 


xo 




xo 


LU 


xo 


zx 


H 


X 


ao 


M 


NO 


< 


i^_j •* •* ii 


i— a 


•»ao 


UJ 


-O' 


X 


•kvO 


1 




•k 


n 


h- 


r- 


CL 


r-OUUQUQGHH^ 


UJ o LU 




X 


— «r^- 


h- 


*-H 00 


<— < Q_ 


H 


<— 4 


•» 


U. 






< 1 <Z<Z<Z< 1<1 II 11 Xuj 


-LCOQC 


II -LU 


H 


II -*LU 




II «kUJ || 2.LUQC 


II 


r — 4 


O 




UJ 


«—*ZUJ LUUU LULU UJ»— ♦ ”0 --J 


H *-4 


OHJ 




x^x 


LJ 


H J 


—•0X0 


4—4 


w 




w 


X 


— ^CLCLCC.OLCLCL mZ 


CL 3 


— 2 


H 




<z 


— X 


II ZOO 




O' 


H 




H 


O ro ms— *— • 


>oo 


O t-H 


UJ 


O H-, 


UJ 


o MO-.M 


1 


o 


O *— ir— « 


X 


a 






LL UJ 


HOh 


OO 


MQh 


CL 


vOQh^nh 


1 


00 


<3 n 


CL 




H 


<Q J-JJ — U J 00Z 


»-»QCQC 


<Z 


1 


<z 




<z 


—2 


: I 




lu M o 


Z 


LU 


X 


lllcL<t <<<<<0030 


o 


oujo 


1 


ouja 




0LU00Z50 1 


o 


c*3Z 


►—4 


cl 


CL 


QCXOOOOOOOOOCO 


LULL 00 


OQCO 


1 


Q0C-O 




Q^OOZU 1 


Q 




1 




X 


1 “*^ 


Q_ *— « ♦— « 




1 




1 






1 






1 

1 




4—1 

1 




00 

1 




1 

1 




1 

1 






1 

1 






1 

1 




1 

1 


1 


1 

1 


o 


1 

1 


o 


1 

1 






1 

1 










1 


CO 


1 


rH 




IN 


1 


o 


o 


► 





ooo ooo 



HUOO oo 



ooo 



ooo 



r-oo 



123 



o oooo ooooooooooooooo ooooo oo ooooo 


00 


OOOOO 






m Qoo , ' 0 ^(NmN 3 -invor^aoa' 0 ^osj mNrm-oi s -coa' 0 ^c\icno 


COvf 


M-^cgrOsT 






CXJ CXj UJ <JU CXJ UJCU LT* L/NLC U' LT*LP» 


J'U’CJUWUUUUUUU 


U 


t-H HHH 






M C\JC\|(\JI\J tNJCNJCM(NC>sjc\|CMCMC\j(\j(\j^\jmnimrnrnrnn^nirnnirr)n^nirnn) 


CM CM 


mn)p)n)m 






O OOOO OOOOOOOOOOOOOOO ooooooo oo ooo 


OO 


ooooo 






Z ZZ 2 Z ZZZZZZZ 2 ZZZZZZZ 


zzzzzzzzzzzz 


Z 2 


xxxxx 






a auuo uuoouaouaauQaua uouoouo uoaou 


uo 


ooooo 






o oooo ooooooooooooooo ooooo ooooooo 


00 


ooooo 








if 




if 














if 




if 














if 




if 














if 




if 














if 




if 














* 




if 








1 






if 




if 








1 






if 




if 








1 






if 




if 








1 






if 




if 








1 






if 




if 








1 






if 




if 








1 






if 




if 








1 






if 




if 








1 






if 




if 








1 






if 




if 








1 






if 


LU 


if 








1 






if 


O 


if 








1 






if 


2 


if 








1 




o 


if 


LU 


if 








1 




o 


if 


a: 


if 








1 




Oh 


if 


LU 


if 








1 




o 


if 


LL 


if 








1 




HO 


if 


LU 


if 








1 




H 


if 


C£ 


if 








1 




o 


if 




if 








\ 




HO 


if 


O 


if 








1 


*-H 


o 


if 


X 


if 








1 


CO 


o 


if 


< 


if 








\ 


in 


o— 


if 




* 








| 






if 


O 


if 








1 


a 




if 


X 


if 


0 






1 


H 


— *r-H 


if 


H- < 


if 


0 






1 




— + 


if 


X 


if 


►— 1 






1 


0 


-HO 


if 


O 


if 


X 






1 


0 


+ — 


if 


IU 


if 




Or 








03 


if 


32 


if 




O 




>• 




—X 


if 


0 


if 


oc 


a 




O 


CM 


00 


if 




if 


•» 


0 


— * 


X 


• 




if 


cc 


if 


0 ^ 


X 


0 — 


0 


O 


— 00 


if 


0 


if 


—I x — — — 0 


•» 


-J oC — • 




UJ 


OC — < 


■it* 


LL 


if 


LUCO — 1 < CO 0 — 1 


r-» 


O 0 LL 


CO 


• 


^ J2 


if 




if 


ohx _j -j 0 a 


II 


ooo 


►— < 


-J 


► <1 >-1 


if 


< 


if 


^>•0 oooo 


X 


x 0 0 


oo 


QC 


0 uj<t 


if 


H 


-If 


O' **oox 0 0 0 X 


•k 


•» X 0 


> 


H 


2 : a: • 


if 


< 


if 


»— CM -X X X X 


— * 


0 - X 


O 


X 


X «H 


if 


O 


if 


CJf 




2 (X 


<1 


a 


•* UJO 


if 




-If 


-^CXC < CO 0 -£ 


•» 


0 3 U. 


X 




0 x • 


if 


H 


mux •'•—3 3 : 3 3 a 




Ct: 0 2 


< 


• 


z •— 


if 


O 


-if a <h 5 -q 0 a 0 a: 




X DC O 




Q 


0 


if 


CL 


ifOZH>a: a: a: or ^ 


03 


- X c£ 


LU 


X 


< 


— if 


z 


if ^ ^ 2 oOZ X X X 


20 


2 - X 


3 


< 


LU -5w 


*"H *-H # 


>— « 


-If 


•*•»•* zs 


a~ 


a 3 ^ 


O 


• 


C* O O'— O 


+ CL if 




if 




— Cri — 


LL 


< 


CM 


— X. 


01 if 


Hi 


if O OOOOO— 0—0 — — 


CMX 


•— ^wO w 


> 


• 


z * *T:ao 


— LU if 


X 


-«- ocnooh^-i <n m nJ- 


m ^ 


cn vt 




O 


QC * — — 'O-'— 


OH if 


H 


if CO l^-in OOOOLUCOLUCOLUCOLUCO^if 


LU CO UJ CO LU 


QC 


LU 


H O ll —0 0 2 - II it 




-If - 




•* 11 » 


1 — »H •‘H 


<1 


• 


Ujul ll ZUJH“)j<ii.u-.UJUJ^ 


LU 


it sU sUsU-Ur'^vO*— <0 »— 'nOh-<sO*— i^U^sL)LU*-HvO*— 




ULi 


Mon-)-* 11 0 l < 3 . 2:0 M '-HOO-* 


H 


-if — 


^wwwwq:wq:wq:wq: 


w w 


ooc—aiwo: 


O 


O 


11 z> h 




►—1 


if 


3333 


ro 


X 3 2 2 


O 


O 


-**-h 0 — o»— •# 


QC 


if LU LULU LULU LU LU LU LU 


ULIHLU 


M LU 111 


X 


*c 


MhJO'Ct^H^HwwQ. 


ww Hh^ 


3 


if hhhhhhoh— IH—IH — 1 




hJh-IHJ 


►—1 




— Z — 1 HDZ 11 




if •-< 


mhhmh Jh 0 *— 1 — J *— • 0 


t-H ►— < 


XO— iO— *o 


CO 




arcxo 11 ?o-)Ou-llujZzoc^ 




•* cc cc cc a: ac cc < a: < cc: < oi < cc Oo; c < cc < a: < 


1 


u. 


hOUQ^QO OQ^^KOOOO-^ 




#S232ii0iU3020 


2Q50U2U30 


1 


»— 1 




if 




if 








1 






if 




if 








1 






if 




if 








1 




t • 


OO 










CO 


1 




0 0 


O-H 










r-H 






CO ON 


HHUOUUUU oooo 00 




HO OOO < 



124 



ooooo 


o 


oo 


uvom<rtr\ 


h- 


00 O' 


«— * f-H r- r» r- 


r** 


r-r- 




m 


mm 


ooooo 


o 


oo 


z :xxxx 


X 


2:2 


UQUUU 


o 


oo 


uuuuu 


o 


00 



I 

I 

I 





1 

1 






1 

1 

1 




1 

1 

1 




t 1 

1 1 

1 t 


l 

l 

l 








1 






1 




1 


O 


1 1 


1 








1 






1 




1 


►— 1 


1 1 


1 








1 

1 






1 




1 


UJ 


1 1 


1 








1 

1 






1 

1 




1 

1 


X 


1 1 
1 1 


1 

1 








1 






1 




1 


UJ 


1 1 


1 








00 






1 




1 


-J 


1 1 


1 








a: 






1 




1 


a. 


1 1 


1 








a 






1 




1 




1 1 


1 




CL 




H* 






1 




1 


5 


1 1 


1 




UJ 




O 






1 




1 


0 


1 I 


1 




►M 




LU 






1 




1 


•«*» 


1 | 


1 




■k 




> 






1 




1 




1 1 


1 




< 










1 




1 


m 


1 1 


1 




3 




a 






1 




1 


X 


t I 


1 




•» 




x 






1 




1 


— 1 


1 1 


1 




O 




< 






1 




1 


tt 


1 1 


1 




» 






*— % 




1 




1 


M 


1 I 


1 




O 




to 


CL 




1 




1 


O 


1 1 


l 




3 




Ui 


UJ 




1 




1 


• 


1 1 


l 




•» 




0 


•—» 




1 




1 




1 1 


1 




< 




0 


•» 




1 




1 


cO 


1 1 


1 




-*3 




<1 


x 




1 




1 




1 1 


1 




<50 




J* 






1 




1 




1 1 


1 




30C 




1 


•» 




1 




1 


O 


1 


1 


<< 


OX 




0 


0 




1 




1 


►—I 


1 i ~~* 


1 


33 


0 ► 




•— • 


H 




1 




1 


UJ 


1 II 


1 


OO 


Zm 




UJ 


•* 




— * | 




1 




1 a 


1 


OO 


►X 






M 




— » | 


x 


1 


-J 


1 XX 


1 


XX 


<3 




UJ 


•» 




•> | 


0 


1 


< 


1 ouj 


1 


•» •* 


3 - 




X 


O 




X | 


0 


1 


UJ 


1 ox 


1 


<< 


OO 






►— * 




X | 


0 


1 


QL 


1 »-iH- 


1 


33 


QL *-H 






UJ 


— * 


a 1 


CO 


1 


w 


1 




OO 


X ► 






•* 


H 


3 1 




1 




1 h-— » 


X 


ora: 


**< 




O 


fNJ 


•* 


•* 1 


X 


1 


V 


1 OO 


>— » 


XX 


X 3 




UJ 


*» 


X 


x 1 


a 


1 


rsj 


1 O • 


CL 


*» * 


CLQ 




to 


O 


(X 


x 1 


*— • 


1 


O 


1 a:o 




•— 1 H— 1 


xo: 




♦— 4 




0 


U | 




1 


m 


1 — * • 


< 


X 3 


LUX 


— % 


3< 


•* 


3 


X 1 


<1 


1 




| — * * X UJ 


X 


3rv| —* 


h“ ► 


0 


< 


<t 


* 


to I 


X 


1 


• t 


| •— OLU X 




► ► LO 


— » -X 


2 : 


xx 


3 


< 


► 1 


CL 


1 


X 


| 0 — 3 • 


x<< 


XX • 


inxo. 


2: 


cl *— « 


O 


3 


< 1 


O 


1 


♦— « 


| —00--* 


-J 33 


33 O 


•3X 


•» 


Oa: 


a:<< 


0 


3< 1 


LL 


1 


CL 


| <0 — X'-* 


OO 


1 


0 UJ 


0 


Zh 


Z33 


X 


03 1 


to 


1 


H- 


I 3*— UJO”D 


U >C£0 


— ■— ► 


h- 


x 


►< 


►OO'-* 


ZSciD- 1 


X 


1 


< 


1 ouj— 0 — 


NZX 


O 


O — 


— * w 


h-x 


<tcLCL'-4 


— DC 


TCC^L 1 


< 


1 


X 


1 X — O O 


►M •» -O O 


OOLH tin *30 


0 


a: 


— xx - 


i\IO 


*X — 1 


oc 


1 




1 X3<TlO*^ 


JMH • • 


00 «o 


•OO-H 


LO CL 


OH 


•> 


XX 


-J -X 1 


K 


1 


X 


1 ►<TX<LU 


< II II O — 1 


UJ UJ 0 — 


O— U!K 


X) 3 


tOX 


LL HHw 


a:to<^x 1 




1 


0 


| huJmXw 


mXJ It II LilJJ II II 


It II 30 


-O 


►< 


cc 11 11 rvi 


ZVUIIU 1 


LL 


1 


►M 


| || CL < O 




OOlO--*-> 


0.111 


vU ULi 


uj_j 


Omi-) II Uio • 


-2UJUJ 1 


O 


1 


h— 


1 3 II II <1 <1 


11 


ZUUh(\)H(MUJ 


^ J> 


h-O- 


*— « —■* 


OtO<-*LO II 30 1 




1 


<X 


1 — XQ-Xlplp ► ► 


M •• •* 


•> •» 


O 


O 


UJHfOOZ || || OLO-*XX | 


X 


1 


X 


1 033LL 


h J3 hh(^jc n | 33 


UJ 


CLUJ 


OOQOw 


hhZT 


ao « | 


•— 1 


1 


CL 


1 QN- — **— 1<0 1 — — 


X33 — 


— 33 


l~ 3 


XX 


3 —i -< ^ 1 — (X or 3 ^ •‘i— h- 1 


0 


1 


CJ 


1 -<<OCO — 1 


0 l OOXXO<<XXXX<<T 


»~i 3 


Ol- 


0 0 : 2:000 ^zz 1 


LU 


1 


LL 


1 MM | 


m | O 03300033 


3 — 1 OO 


c£ < 


0 


< oooazz 


< 0-00 1 


CO 


1 


to 


1 OUJlil 1 LL 


1 






-20 


1 LL 


UQQiU^UQNUU 1 




1 


X 


I oa:^M 1 m 


1 








1 O 






1 




1 


< 


1 1 


1 








1 






1 




1 


oc 


1 1 


1 








1 




m 


LTW t | 




1 


h- 


1 1 


1 


h- 






1 




HD 


com 1 




1 




1 1 


1 


ON 












3.-J 1 




1 




1 




H 





OOoO ooooooooo o o o 



125 



0 


000 


00000 


00 


OOO 


0 


-Hcxicn 


<j- m so r^-co 


00 


c^o>a> 


UJ 


UJ UJUJ 


uj aj uj ujuj 


UJ U' 


u \u \U\ 


cm 


cnnim 


nimrMmo) 


u)m 


CMn>m 


0 


00 0 


00000 


OO 


OOO 


2 


zzz 


Z2ZZZ 


22 


zzz 


0 


uoo 


uauuu 


OO 


OOO 


0 


000 


00000 


OO 


000 



o 



U-» 



O 3 IIOOOOO 0*100 

ZQUJUJUtoO^>^>Ow-)0 --ZZ 

• ^h— H-H- O O c\j> — ~«ax»~<»— • 

^ — l_l» * — * . — » » — » ( rvi i — ic\jarj^^Q.I — H- 

2 1 CLQCCC-l >-J-K- O 2*-« — 222 

i .^2i<ouj<^<a iiuioniuoo 



a 

o 



OO^o# UQ^OQOHUU 



LL 

O 

l-U 

o 

< 

o 

a 

2: 



1 






1 a 


1 




1 




1 






1 uj 


1 


CC 


1 




1 






1 X 


1 


— LU 


1 




1 






1 


1 


oC •— 1 


1 




1 






1 a 


1 


LU - 


1 




1 






1 UJ 


1 


*-*< 


1 




1 






1 H 


1 


-3: 


1 




1 






1 or 


1 


c •» 


1 




1 






1 UJ 


1 


30 






1 


m 




1 -2 


1 


•* ♦» 






1 


it 


it 


1 *”H 


1 


00 






1 


it 


it 


1 


1 




LU 




1 


it 


it 


1 UJ 


1 


O - 


O 




1 




it 


1 Cd 


I • 


^C 


O 




1 


-J 




1 


1 it 




c 




1 




LU 


1 > 


1 it 


— CO 


> 




1 


z 


O 


1 <1 


1 it 


C iQC 






1 a. 


0 


O 


1 5! 


1 


— O ox 


0 




1 2 


M 


< 


1 


1 < — 


CO X - 


►— f 




1 0 


H- 


> 


1 LU 


1 < 


00 2-« — 


LU 




1 0 


< 




1 2 


1 LL.O 


02 -I— O 






1 •— 




oO 


1 ' 


1 OO^ 


O - 


X 




1 •> 


a 


QQ 


1 H- 


1 o< 


Z< NO<cQQ 


LU 




1 • 


0 — 


< 


1 O 


I a :2o 


-3 O JU 


O 




1 


LLC 


— * 


1 O 


l 0 -0 


CO •‘O 0 0 2 


a 




1 * 




>a> 


1 CC 


1 HCo 


2qT O *— <00 °* 


2 




1 it 


Z3 


luouj 


1 


1 0X2 


02 -h *>zza 


O 




1 * 


<0 


bZ 2 


1 H- • 


1 UJO * 


a: - -< - -3 


O 




1 


OCX 


IN 


1 2m 


I >ac 


2 (XJ C 2<CQO 






1 • 


H* •* 


cac 


1 ujo 


1 ZZi 


-X 3022^ 






1 a 


< 


jS CD 3 


1 2 : 


1 u.i -a 


~a 02002 


X 




1 2 : 


X3 


Qaiu 


1 2X 


1 Ocx a: 


H-2 axaa * 


c 


2 


1 0 


OO 


a: oar ~* 


| O — 


1 *-»CX2 


NUJ 2 -223 




LU 


0 


car 


202 2 


1 •—'or 


1 1112 - 


OH- -CXJ - -O 


2 


X 


Q.*—» 


2: 


•>C •* C •* 


1 </>H- 


1 ujo: 


- -MXXCcO - 


a 


H 


O 


it •*< 


O OC— 


1 uoc 


1 it H-a<< 


hNOO-Q. - -X 


0 




Ott 


# xz 


»-<>-— 0-) 02 x 1 


1 <Z 


1 # CZJ 


mo zrxxz 


a 


— N 


-J* • 


- -JO 


uj cq uj 0 1 —a: •» — 


1 


1 ft nludo 


O >LULUCcQO 


00 


•—4 




cc 


c£ •— « O ^ 2 2 2 M 


1 ala' 


1 - Oharu 


—— 0 X 1 — 1 — - — ' — — 


2 


• 


O'"* — 


2 


— O'— 2 + — -X II 


l 00 


1 —— 2 : 2:00 


u— — 


C 


O 


Ott # * LU - 


2 -<o^— * 


1 H-h- 


| -N<UJ *‘N|I 


DDhUO»> 


a 


LU 


•* •* * 4— «— 4 


a*-»a — «h 11 m 0 


1 00 


| ■#■ UJhHHOH 


oo^^^ooo 


H- 


• 


U. O'L 


tl 


11 aiux 11 - 


1 Lum 


1 »a»-i i| 11 11 11 


UJUJ UhHXXX 




a 


LJ — 


- — (X O JQL (X LX. —)LL. II *;-n*.UJLU 


1 >> 




UJJJUULPJJJ 


1 


2 



22 

ILJUI 

OO 



UJUJ 

>UJ 

2X 



— Q — £ 0-«Oaa 

o cxjcxi - * 2t 
LU (XJ (XJ *— •»— 1 »-i ' 

1 — 1 i - — ' •— f — ^- 1 — n — — 1 — »— 1 — 1 

— 1 — I — J — IC2? — I — » — I — I — » 

a'<<OCNNC<<tLL<<<<I< 
SUUQOOOOOOJUOUUU 



O l 


0 


00 


(XJ 


O' 1 


0 


1 — 4 (XJ 


(XJ 


1 


(XI 


NM 


<\J 



O 

o 



OO 



000000 



OOoO 



126 



I 



I 





1 

1 


cC 

UJ 


1 

1 






1 

1 




1 

1 




1 

1 




1 

1 


1 

1 


~<\l 


1 




1 


XX 




1 




1 




1 




1 


1 


XX 


1 


•> 


1 


a.3: 




1 




1 




1 




1 


1 


CLO- — 


1 


< 


I 


20 




1 




1 




1 




1 


1 


-*22X 


1 


->3-* 


1 


— UJ - 




1 


— * 


I 




1 


— — «* 


1 


1 • 


cmujuk 


1 


CM -X 


1 


c\lh-< 




1 


— X 


1 


X 


1 


XX 


1 


1 


Xhh - 


1 


XOCO 


1 


X -3 




1 


x< 


1 


CO 


1 


CL JC 


1 


1 ■# 


CL - -< 


1 


CL — * 


1 


cl<o 




1 


CL - 


1 


•> 


1 


20 


1 


1 •* # # 




1 


2 00 


1 


23: QC 




1 


2 < 


1 


CO 


1 


LU - 


1 


1 # # # 


UJ220 


1 


UJ «— *_J 


1 


U'QZ 




1 


UJ3 


1 


O 


1 


h~< 


1 


1 * # 


hOOa: 


1 


h— — O 


1 


h--cL - 




t 


H-O 


1 


a 


1 


-3 


1 


1 X 


-QLCL2 


1 


-<o 


1 


*20 






-CL 


1 


o 


1 


<o 




1 X X 




1 


<32 


1 


<x -o 






<2 


1 


2 


1 


3CL 




1 CL ♦— 1 


“5 - -< 


1 


-SO - 


1 


-SC JO 




UJ 


3 - 


1 


•. 


1 


02 


o 


1 h— CL CL 


o«j 


1 


-*0 QC < 


1 


Ooo 




o 


o< 


1 


— < 


1 


QC - 


o 


1 < 1 h- 


q:jjd 


1 


xcl23 


1 


CLQ2 




o 


CC -i 


1 


CQ3 


1 


20 


o 


1 2— <0 <— 


2000 


1 


jz -o 


1 


20 - 




< 


20 


1 


OO 


1 


-O 


CO 


1 < 20 2 CD 


•*002 


1 


o -—a: 


1 


-2 0 




> 


-O 


1 


CDCC 


1 


oo 




I <o o O 


<22 - 


1 


0<02 


1 


< -3 






<2 


1 


02 


1 


oo 


LU 


0*00- *CQO — 


_J - -< 


1 


20^1 - 


1 


ooo 




o 


o - 


1 


2 - 


1 


02 


LO 


00O< 20 OCO 


0<< z 


1 


-O -< 


1 


O30C 






o< 


1 


-< 


1 


O - 


o 


-OJZOQ -0020 


03 30 


1 


CG003 


1 


002 




LU 


03 


1 


C Q2 


1 


2Q 


LU 


CO 2 -o LUO OLU **o 


zaoo: 


1 


2ZHQ 


1 


2cL - 






20 


1 


30 


1 


-3 




-a: <0 2 30 2Xo 


-rLcL 2 


1 


O - **cL 


1 


-2X 




o 


-CL 


1 


OcL 


1 


oo 


2 


<032cL02lL32 


<22 - 




aL<<2 




< -CL 




< 


<2 


1 


lL2 


1 


3cL 


LU 


LL O -O CL -OO - 


.2 - -og 


3 


Z3ti - 


X 


3: CM 2 




LU 


5 - 


1 


2 - 


1 


02 


X 


OL/)CL<LL20LLCLX 


OMXX 


CD 


-aox 


2 


OXUJ 




QC 


OX 


1 


-X 


1 


CL - 


I'- 


DZZ3C0 -3U023 


QC X CL CL 




XCLCLO. 


O 


OCCL h~ 




w 


CL CL 




XCL 


1 


2 X 




2 < -ozccaz -o 


2CL2 2 


2 


CL 2 22 




22 - 






22 


X 


0.2 




-CL 


LL 


qcclclcl<o-cl<clq: 


•*2 LULU 


OC 


2 - -LU 


2 


•.LUCNJ 




X 


•.LU 


x 


2 LU 


X 


N2Z 




Ol— CL 2 CL 2 2 CL a. 2 


XLMh“H 


o 


aJxrM»— 


uL 


XhX 




< 


NJI — 




UIH- 


3 


OLU 




LL X — f— LU **h-2 - 


-Jh- - - 


LL 


h- ox - 


o 


O -CL 






O ~ 


2 


h- - 


o 


-h- 


LL 


LO-^ UJ CL \—cC UJ CL 


•» i *— i 


O0 


- **CL ►— i 


u. 


*X2 




2 


•»*— « 


CC 


t 




X - 


o 


2 ^- | — CL -CL & h— CL 




2 


XM20 


uo 


N"?UJ 




QC 


XM 


a 


XN 


2 


3M 




<tH -2^X2^h -21 


o<oo 


< 


CQOLUt-M 


2 


OOh 




o 


<o 


LL 


COO 


cc 


OO 


o 


CL XLU# 3LU-M* XUJ 




cl 


w w| 


< 


W-ww 




U. 




U1 


—• — * 


o 


WW 


2 


- <>— - Oh-- XH- 




h- 




m 






oO 




2 




U- 




LU 




JJJJ 




DJOJ 


h- 


OOh— 




^ — 


OO 


< 


0-1 


t/) 


OJ- 


1 


h— <UJ— <UJ-»<U1 


I5III 


i 


02 <— »2 




222 




< 


22 


CL 


02 


2 


22 


1 


o* UJ»— LUj—ih UJH- 


H— t— 


i 


U lh-h-h- 


1 


f — h— H- 




□1 


H-h- 


)— 


LUt— 


< 


h-h- 


1 


Q- —CC>^ *-CL ♦ —CL »— • 


<L<1<<* 


l 


-J<1 o<i 


1 


<<< 




1 — 


<A< 




o<a 


iX 


<< 


1 


h—MJOLLOOQC'OOCL 




i 


Q-2LU2 


1 


222 






22 


1 


CL2 


h— 


22 


1 


j-q.j:w(lj:wq. jc 


oooo 


i 


OO-JO 


1 


OOO 




I 


oo 


1 


OO 


1 


OOLL I 


o o o o 




i 




1 






1 




1 




1 




— 1 


I LU LU LU 


1 JO I 


l 


OOOO 


1 


OOO 


LU 


1 


OO 


1 


oo 


1 


OO 


1 


| | OI—OOJ IO 


oooo 


i 


-JO-JO 


1 


ooo 


oo 


1 


OO 


1 


oo 


1 


OOO | 


1 *— •OO’— ‘00*~*00 


<<<< 


| 


<<<< 


1 


«< 


o 


1 


<< 


1 


<< 


1 


<< 


2 1 


1 CL<<CL<<CL<< 


oooo 


1 


OOOO 


1 


ooo 


UJ 


1 


oo 


1 


oo 


1 


OO UJ 1 


1 300300300 



I I I I I I I 

I I I I I I I 

I I till 

1*1 I 

o ouu ooo o uuu ooo ooo OOO 



127 



END CF TRANSFORMATION BLOCK 



oooooooooo o 
HNrn'4*inor^coo>o o 

SJ v(j \U U 'Ll r— vU 

CNJtNjrvjrsjCNUM C\J^ m 

ooouoouooo o 

X 

u 



auQuuuuauu 
oooooooooo o 



o 

O' 

ay 

<3 

O 

2 ! 

O 

o 



00000000000000 
r* coc^o-HrMm-r^or-co^o 

r-l r-« ~ 4 ( \JlM l \4 l \| l\J l \| t \g < \J t \j ^ \j ( t j 

333333330)330133 

OOO OOO OOOO OO OO 

ouuouuuaauuuua 
OOO 000 0000 00 00 











1 

1 1 






1 

I 




1 

1 




0 

•» 


1 








1 1 






1 




1 




CQ 


1 








1 1 






1 




1 




3 


1 






1 


1 1 






1 




1 




-> 


1 


1 




1 


1 1 




Od 


1 




1 




•» 


1 


1 




1 


1 1 




LU 


1 




1 




CQ 


1 


1 




1 


1 1 




rvi 


1 




1 




O 


1 


1 




1 


1 1 




hH 


1 




1 




> 


1 


1 




1 


1 1 




X 


1 




1 




•* 


1 


1 




1 


1 1 




►— 1 


1 




1 




X 


1 


1 




1 


1 1 




H- 


1 




1 






1 


1 




1 


1 1 




a. 


1 




1 




X 


1 


1 




1 


1 1 




a 


1 




1 




0 


1 


1 




1 


1 1 






1 




1 




0 


1 


1 




1 


1 1 




UJ 


1 




1 




X 


1 


1 




1 


1 1 




X 


1 


00 


1 




•» 


0 


1 




1 


1 1 




t- 


1 


•* 


1 








1 




1 


1 1 






1 


a> 


1 


-J 


Qbd 


0 


1 




1 


1 1 




00 


1 




1 


<r 


X3 


X 


1 




1 


1 * 1 




CL 


1 




1 


3 


**— 1 


t— c 


1 




1 LU 


1 # 1 




X 


1 


00 


1 


X 


Qod 


H 


1 




1 X 


1 * 1 


* 


3 


1 


LU 


1 


< 


<X 


od 


1 




1 M 


1 tt 1 


-* 


3 


1 


CQ 


1 


X 


Qd •* 


< 


1 




1 H~ 


1 1 


# 




1 




1 




o^d 


H* 


1 




1 3 


1 -* oc \ 


LU 


O 


1 


LU 


1 


LO 




00 


1 




1 a 


1 0 1 


ex: 


LO 


1 


X 


I 


a 


*+— « 




1 




1 oc 


1 X | 


LU 


< 


1 


h- 


1 


< 


H- * 


LU 


1 




l 


— - 0 1 


X 


CL 


1 




1 






X 


1 




1 0 


r-* 0~* I 






1 


Qd 


1 


Ul 




H- 


1 




1 LU 


U OC I OC LU I 


*: 


H* 


1 


O 


1 


UJ 


CdQd 




1 




1 INJ 


LL IU-JOO | 


0 


(3 


1 


UL 


1 


Ul 


CL X 


X 


1 




1 ex: 


x x • ood 1 


0 


ex: 


1 




1 




t— 4 * 


O 


1 




1 3 


HH *—*•—< OU.1 | 


0 




1 


LU 


1 


0 




oc 


1 


— * 


l s. 


* oc | 


CQ 


LL 


1 


> 


1 


< 


OX 


LL 


1 


oc 


1 


~i H- ON *X \ 






1 


-J 


1 


-J 


LU ."5 




1 


0 


1 0 


m <c'od*-‘ 1 


~rr 


LU 


1 


O 


1 


LL 


X * 


>< 


1 


0 


1 X 


x — 1 


O 


X 


1 


OO 


1 




a< 




1 


0 


1 c 


oh oahec 1 


t — « 


H- 


1 




1 


LU 


^0 


UJ 


0 1 


2 : 


1 


x - x +cc oc: 1 


H- 




1 


a 


1 


-J 


-X 


0 


3 1 


•* 


1 UJ 


OH- OoCZOO 1 


< 


O 


O I 


H- 


1 


< 


H- * 


CO 


1 


oc 


1 X 


OUJ 5 -cCcC 1 






3 1 




1 


O 


CL< 


< 


O 1 


JZ 


I ►— * 


qlQ h-o oc X 1 


►— t 


2 


-H | 


QC 


1 


00 


ao: 




H- 1 


0 


h— 


• * OQC.Q-LU » I 


X 


< 


1 


LU 


1 




»— X 


Qd 


1 


OC 


NG 


cccc a z^coccc 1 




X 


O | 




1 


LU 


«b •» 


< 


O 1 


X 


XCJ 


-jo a: *uj <3 1 


H- 


1— 


1— 1 


►— t 


1 


X 




> 


O I 


•» 


hq: 


oa cx:h- 3 o 1 


CL 




1 


X 


1 


1— 


<3X 




— | 


_s 


1 


00 ujcl »uce 1 


U 


ex: 


a 1 


*-4 


1 




Qd * 


X 


— ^ 1 


od 


* ^ 


xx ex: xx 00 2 1 




LU 


0 1 


H- 


1 


00 


H-X 


O 


O - 1 


•» 


* 0 


- - <iuar*|f • 1 


Qd 


H- 


1 


CL 


1 


H- 


00 Q 


•— * 


• O I 


X 


oc < 


oooca: 3 HLO^a: 1 


LU 


< 


1 


O 


1 


LU 


f—i •* 


00 


OOO — 1 


3 : 


a. 


-h — 0 *:2 -h-cl 1 


H 


LU 


0 1 




1 


00 


*0 


UJ 


3Jlii x 1 


ex: 


LU 2: 


- **003 co^odtt x 1 


X 


Qd 


• 1 


00 


1 




o~ 


O 


00 • *2 1 


w 


H*-h 


^3:cccc * ex:—* xi 1 


UJ 


O 


1 - 1 


3 


1 


0 


LL * 




CCU-' 0 1 




<J 


occczz*-* UJ— 1 — 1 






0 1 


<3. 


1 


0 


ZH 


LU 


1 


> 


O 


•» «>W f— 1 




H- 


• 1 




1 


<c 


HO 


X 


* ~ **d-J 1 


0 


30 


Omhhj: oluo-R'Uj 1 


# 


LU 


-j 1 


LU 


1 


0 


— X 


H- 


-H^3 *<X | 


X 


02: 


OCQ II II CC CL^VIX *h- 1 


* 


L/ ) 


ad 1 


X 


1 




•* 




II II —■ 3 UJLULUUJ | 


0 




00>-»3 11 LU ex _J'U *— • | 


■* 




h- 1 


H- 


1 


f— 


OOO 


LU 


| 


CL 


<5.00 


cl cl ^ouoa.— oc 1 




-J 


2 | 




1 


00 C\J OO 


H> 


Oi 11 2 : 2 : 2 : 1 


O 


03 


on^HHMZOiO 3: 1 


* 


DC 


O 1 


-J 


1 


Qd 


1 <M 


< 


OOhQ^h m m 1 






33 **-h 1 LU 1 


# 


1- 


X 1 


-J 


1 


t — t 


II 


O 


1 


-J 


1 


JJHHOh 1 OOH--J 1 




2: 


— 1 


< 


1 


LLO OO 


O 


11 :xxxx 1 


O 


1 


-JO —X | OOH-iO I 




O 


1 


O 


1 


l 


xo * 


X 


OOiL'^wOOG | 


< 


1 


<<0020 | «a:< 1 




^d 


UL 1 




1 


O02<3 


a: 


QQh^xUOO 1 


O 


1 


ooaoaio I 00 50 1 






M | 




1 


hmOQ] 


a 


1 




1 


1 1 




LL 


1 




1 


X 


»— H 


ll 


1 




1 


1 1 




♦— H 


1 




1 


H- 




1 


1 




1 


— » 1 1 






1 




1 






1 








3 1 






1 




1 








000 1 






r— < | 

















000000 



rO'tLHO OOOO 



OOOOOO OOOOOOO 



128 



o 



ooooooooooooooooooooo oooooooo 

^or<0 CO O' 'HOC'grO ICO Ob-COC'O-H CMm nT lAsO 0^0 -h -H 

niiMiM(^M\in»oi(MnnN4>r '0'NTiN>r'r'J''T'Tma» lmu*u munutmiM 

rncnnqqrntTvqrncqcocqcqcnrqcqrncncqcncqrn cqiqrqcqnHMrqcn rq 

ooooooooooooooooooooo oooooooo o 

22ZZZ22ZZ2ZZZZ22Z2Z2Z ZZZ2ZZZ2 2T 

uuauauaauuuuauGuuoaua uauauaou o 

OOOOOOOOOOOOOOOOOOOOO OOOOOOOO o 



1 






a i 




1 






1 




1 






1 






•* i 




1 


1 




1 




1 






1 






CD | 




1 


1 




1 




1 






1 






ID | 




1 


1 




1 


CM 


1 




1 


1 






> | 




1 


1 




1 




1 




1 


1 






9* 




1 


1 




1 


pH 


1 




1 


1 






QJ | 




1 


1 




1 


1 


1 




1 


1 






U | 




1 


1 




1 




1 




1 


1 






> | 




1 


1 




1 


if 


1 




1 


1 






•» | 




1 


1 




1 


if 


1 




1 


1 






X 1 




1 


1 




1 


CL 


1 




1 


1 






•* | 


X 


1 


1 




1 




1 




1 


1 






X 1 


o 


1 


1 




1 


OO 


1 




1 


1 






a i 


*— ♦ 


1 


1 




1 


w 


1 




1 


1 






o I 


1 - 


1 


X 




1 


o 


1 




1 


1 






z i 


< 


1 


»— ♦ 




1 


if 


1 




1 


1 






•» i 


h- 


1 


cL 




1 




1 




CD 


1 






>— i 


3 


1 


b- 




1 


o 


1 




* 


1 






□ ^ I 


CL 


1 


< 




1 


* 


1 






1 






^2 1 


X 


1 


5 : 




1 




1 




"H 


1 






♦ | 


O 


1 






1 


b- 


1 




1 


1 






Qoc 1 


o 


1 


X 




1 


if Z 


1 






1 






< 2 : i 




1 


LU 




1 


if Q 


1 




* 


1 






G £ - I 


—1 


1 


H" 




1 


— t 


1 




if 


1 


00 




O* 1 


< 


1 


00 




1 


O 0 h- 


1 




— X 


1 


00 




•-•2 1 


3 


1 


>- 




1 


1 O 


1 




<a 


1 


UJ 




•*»“H j 


b- 


1 


00 




1 




1 




i 


1 


o 




b- - | 


O 


1 






1 


03 


1 




• b- 


1 


O 




1 


< 


1 


o 




1 


— LL 


1 




O 0 < 


1 


CC 




-•2 1 




1 


LU 




1 


CM 


1 




w O 


1 


CL 




oca i 


LU 


1 


CC 




1 


^LU 


1 




OU 


1 






OLZ 1 


X 


1 


►— i 




1 


o> 


1 




II O 








►-H •* | 


b- 


1 


oo 




1 


1 -• 


1 




o 


o 


o 




j 




1 


IU 




1 


wh 


1 




LU 




1 — * 




QX 1 


00 


1 


o 




1 


* o 


1 




LUO 


*— « 


b- 




id 3 : i 


LU 


1 






1 


if LU 


1 




CC O 


O 


< 




z: - i 


O 


1 


UJ 




— 1 


QCO 


1 




UJ CL 


< 


rM 




o< 1 


Q 


1 


X 




OIL 1 


+ aQ 


1 




X 


O 


♦-H 




H(J | 




1 


b- 


O 


-J — 1 1 


-*a 


1 




20 


oo 


3 : 




* 2 : i 


X 


1 




<r 


au i 


W 


1 




UJ 








b- - 1 


< 


1 


X 


pH 


OO 1 


00 


1 




— a: 


o 


b- 




CL< | 


a: 


1 


h- 




2 ZX 1 


1 — • 


1 




oo 


x 


CL 




OalOO 1 


o 


1 




o 


•> »- j 


UJOC 


1 




woo 




O 




■-Hizaor- 1 


a 


1 


2 


b~ 


OU. | 


QLU 


1 




OLD 


oo 






•» -<r nt i 


CL 


1 






22 | 


IM 


1 




a 


UJ 


UJ 




i— <t i 


CL 


1 


X 


O 


OO | 


UlH-l 


1 




LU 




X 




<xoa i 




1 




CD 


CLCC 1 


XX 


1 




f— LU 


o 


1 — 




DC *b-|- | 


LU 


1 


QC 




ZZ 1 


1 *' * 


1 




3 X 


CL 






1 — X 1 


X 


1 


o 




*k •» | 


b— 


1 




0 . 1 - 


21 


00 




OOQOQ | 


1 - 


1 


1 “ 


o 


OIL 1 


LUCL 


1 






»— < 


b- 




»— t »*00 1 




1 oar 


o 


• 


•» •» | 


b-a 


1 




OoO 




CL 




•*o | 


X 


1 20 


UJ 


or 


XX 1 


3 


1 




Om 




< 




o-»— - 1 


o 


1 OO 


> 


UJ 


OLL 1 


CLUJ 


1 


3 






b- 




X **0 <— H 1 




1 CL O 






w PW 1 


XX 


1 


X 


phOO 


O 


oo 




Zb- • • 1 


z: 


1 zz 


X 


— 


1 


UH 


1 


X 


t— 


m 






«OOh 1 


a 


1 *» •» 


o 


x—xx 


>> 1 


o 


1 


•» 


ZLU 


O l| 


— 1 




— z:luo i 


t—t 


| ,-H pH 






oo 1 


X 


1 


pH 


<X CC 


• 


O 




O • « 1 


b- 


1 II II 


oo 


^3 •» till 


XX 1 


O— t 


1 


II 


-Jill 


0—0 


<1 




ooouu 1 


CL 


1 ox 


LU 


'-'-O — — LULULU 


_J 1 1 


1 — 


1 


1— 1 


ax 


1 <M II 


o 




□ ULUL | 


O 


1 


OJ 


2— 3^333 


CL CL 1 


o 


1 




2 


II » — 






<hZZ I 


CL 


1 oo 




02 - -ZZZ 


oo 1 


CLUJ 


1 


OO 


1 




00 




••KHM I 




1 som 


UJ 




1 


Ooo 


1 


o 


1 


rq3:w 


>-H 




O Oww | 


00 


j f-H *-H 


b- 


'’-**—* ^ w b— f~* b— 


OO | 


03 


1 


<M 


1 




X 




—1 •* 1 


»— i 


1 


<1 


II xxzzz 


OO 1 


o 


1 




1 


* 1 -2 


1 — 




<OlUL | 


X 


1 OO 


3 


H^3.iOQO 


<< 1 


UJ 


1 


O 


1 


-2 1 -* 






OCQt-iM | 


b* 


1 QQ 


O 


H^OOUUO 


oo 1 


XaD 


1 


O 


1 


1 






H j 




1 


LU 




1 


i— i 


1 




1 


1 






1 




1 


1 




1 


<u 


1 




1 


1 






1 




1 


| 




1 


Xb- 


1 




1 








O , 1 




1 O 


1 


ooo 


1 




1 







<m m i vftfvo 

O OOOOO OOOOO^ OOO — < O H(J ooooooo oooo 



I 

I 

ooo 



129 



CALL PL ANTI (GPS, AX , 8X,CX,MU(I ) ,NROKA,NCOLA, NRGhB, NCOLB , NROWC 
1NCCLC) 

PLANT1 COMPUTE G ( — S ) WHERE G = C( S I-A ) ** ( -1 ) *B 

WHERE S IS THE DESIRED POLE LOCATION 





oooo 














oo 




OO 




oooo 






V^- «4* n} - 














O H 




CO 2* 




X)3Q-< 






M(\j 














oo 




— — 




— <•— « \J 






-rvTsrvr 














u^u> 




mm 




mau ujo 






oooo 














OO 




oo 




oooo 






zzz/ 














2 2T 




orz: 










aouu 














CJCJ 




oo 




ouuu 






oooo 














OO 




oo 




OOOO 

1 




•» 

o 










1 

1 




i 

l 

i 


1 

I a 




1 

1 

1 

1 




1 

1 

1 

1 


1 

l 


* 










1 




i 


1 UJ 




1 




1 


1 


a 














i 


1 X 




1 




1 


i 


a: 










2 : 




i 


1 — 




1 




1 


l 


27 










o 




i 


i n 




1 




1 


1 


•» 










*— < 




i 


1 o 




1 




1 


l 


an 










h- 




i 


1 U< 




1 




1 


l 


_j 










O 




I 


1 X 




1 




1 


l 


O 










O 




i 


1 


XICNJO 


1 




1 


l 


O 










o 




i 


X UJ 


■* ir Hr 


1 




1 


l 


Z T 










LI- 




i 


L U — 1 


* * # 


1 




— H 


1 


•» 














i 


MO 




1 




to 


l 


an 










LLI 




I 


— CL 


— 


1 




1' 


1 




— fc 








> 




l 






1 




a 


i 


o 


o a: 








t — * 




i 


— X 


~h c\J cq 


1 




n 


1 


or 


to-^ o 








H- 




l 


H- — 


www 


1 




X 


i 


3 


XX O 








O 




I 


X< 




1 






l 


•* 


to O—O 








LU 




l 


ox 


XXX 


1 




X 


l 


< 


CLOX x 








3 




i 




LULUU! 


1 




LU 


l 


O 


o **o •* 








an 


— o 


i 


U-X 


j— i — i — 


1 




-J 


i 


n 


-xox 




— 




o 


2 •— < 


i 


OlU 


ajiuu i 


1 




X 


LU 


o 


UlOO-J 




— - 






Cl - 


i 


o 


aoo 


1 






X 


3 


-JO -o 




— * 




Ul 




i 


<x 


WWW 


1 




o 


l!J 


•* 


ooxo 




N 




X 


J2x 


i 


*-*x 


OO 3 


1 




o 


X 


<r 


oezoz: 








h~ 


-a 


i 


xo 


<<< 


UJ 








c 


x -o - 




H 






O Jt 


i 


UJO 


XXX 


X 




X 


o 


a 


-ooa: 








—r 


2 - 


i 


h* 


( — 1 ♦— 1 * — • 


LU 




CJ 


UJ 


X 


< jZTj 




1 




»— l 


ny — * 


i 


•— TT 


<C<t < r-4 


X 




U- 


o 


ZT 


jo -a 










OUJ— ■ 


I 


XO 


w w w ^ 








3 


•» 


JO<U 




tt 




o 


O.ON 


i 


0<T 


+ > -h X 


o 




M 


3 


— 


OZ 




x 




u; 


«. w 


i 


LU 


r\ic\j(Aj - 


UJ 




| 


O 


»— i 


2 : -a ~ 




— 




—* to 


h-h-K- 


i 


h~ 


*i\ K — 


1 — 




II 


X 


w 


-OO CT 




— * 




X 3 


> >U J 


i 


to - 


K K -H 


X 






h- 


3 


< *2 ^ 




to 




•* 


CLLXQ 




OH 


% w 


LU 


CM 


33 


r 


x 


>o *o 




— ■ 




3 LU 


» — « ►— 1 '— *• 




o to 


— — — X 


to 


) — 


CUw 


»— 4 


i 


Oa'<^ 




3 




— iGLU 


•» •* — J 


CL 


o 


fX 


z: 


to 


30 




*+ 


XO^Z* 




— 




X X 


Cl lC <T 


O 


UJo 


-H'NJCOUJ 


►— * 


a 


-O 


UJ 




2 ►O •* 




o 




(3 Q — 


O-JUI 


a 


X 


www | 




o 


m »-H (~Q 


an 


o 


~GJXX 








O 1 — 1 — 


□ UlL 


-j 


h-Ui 


Zli-ZUJ 


o 


a- 


w O 






D 




h- 




CO 3 


oo ^ 




X 


XXX O 


►— 1 




2* 


z: 


X 


OCi **o 




■fr 




(X h~u 


z: 2 -5(“ 


u 


2ZU — 


LUUJLU - 




h- 


27 


<T 


d) 


f>0 tO X O 








+ xx 


•■ •> * 


CO 


»— • CNJ 


1 — 1 — 1 


3 


to 


do 


O 


•* 


£0.0 - 




— * 




— < 


ooo 




Qh 


LL'LJUIH 


o 


o 


— H* 






3COOX 




'/) 




x zx 


— 4 — 4 


li- 


UJUJtO 


nn o to 


►— i 


o 


h— 


to 


< 


•* - **C5 




1 




•* mJ 


- «H<- 


en 


U0XO 


WWW CJ 


!- 


+ 


<o 


h- 


•» 


tO X tO tO 




w 




X of 3 X<r 


xo — 




3 — 0 


OO _J o 


o 


M 


330 


zr 


VI 


o. c^:_w 




O 




-£-J — 2 CL 


OQH 


o 


O 


<2 •* 


27 


1 — 


COO 


»— < 


x 


OO oar 




— 




OD — LUZ 


Cl 


2T 


OOQ 


LUUJLU • 


3 


to 


aiu — 


< 


3 


W' ' w 








xo x i — 


*— '►-<! — 


u I 


f-U2 


XX X X 


IX 


3 


- o 


X 


— 






<\l 




Z^OUA LU J 


X XUJ 


1 


x <r 


WW W » 




O 


h- *c\j 


i — 


— < 


JJh _J 




X, 




• - • • U Q 


Q 


1 


to 


WW W 


LU 


+ 


— :ror -o 


to 


h- 






*— 4 




o -h O ^ * D 


O II 


1 


UIU— 


— 


> 


—4 


^lUHQ 


3 


z 


H- h— 




• — - 




ii ii m ii x 02 * 


on — 


1 


1 — ' 1 1 — 




— 


1 — * 


-<t • II — CMO 




<1<1-X<3 




1 




X MJJUIO — 


LUUJ»— * 


UJ 1 


— LJ t > 


•» 


1 — 


to 


nUX X 3 2.UJ I 


o 




i ^ 




•tf 




xo nJ J— 


t JO — 


O I 


h-xo 




o 


o 


wh-U II 3 II 




a. 


-'OOOO 








~4 0 - ~w 2 : 2*00 


O 


2 1 


o 


II II 11 — 


LU 


o 


tOOO — 3 — 


X 


i 


o 




X 




— w — 1 ~5 — 5 of — «*— • 


X 


| 


<0 — 


—c\j n uj 


3 




u'zr^-n 3 — niu 


o 


— 1 — 1 — 1 -JO 




+ 




I 


-j jjj 


1— 1 


3<l- 


H-r-h f— 


□J 


II 


K“ O w vj* w^—< w 


X 


o' J o x o o 




> — t 




2 1 


_J Jh- 


2* 1 


L JIULT) 


to to tO — 


CU 


3 


— O C327C0h- 


< o <f <t <X <X 




— 




OCXX-HfO 1 


<r <r(u 


•3 1 


XO 


one- x 




an 


X LLOOOC 


-a 


ozuaoo 








QOXXXOO 1 


OoO 


O 1 


X. o 


OOO £ 


LU 


o 


jC | mQmOh | 




— ♦ 








1 




1 


x x 




X 




1 


i 












1 




1 


ooo 




i — 




1 


i 












cO — « I 




1 


X -2 3C 




i 




1 


i 
















n i 


1 Li H— 




i 




1 o 


i 












r> rjs 




o 


1 




i 




an 





O'JUUuJOU JUOtJU UUO <r O 



UU' J 



ooo 



130 



G( U = SLFTU REAL (DETERMI1) )!**«' + l A1MAG(DETERM(1) U.05 CGN0E24 

G( 2 J = RE4LCLETERM1 in-C.9 CCN052h 

G( 3 )= X ( 1 ) -X ( 2 ) 



oo 



ooooooo 


o 


O 


o 


or- coo 


ro 


>— 4 


OsJ 


tNiNiviiMCMu i 


(M 


a t 


u \ 


snuunuuninin 


UT 




in 


ooooooo 


o 


o 


O 


zrx-xxxxx 


X 


X 


X 


QUUUUUU 


o 


a 


o 


00000(00 


o 


o 


o 







1 

i 




1 if 
1 - w 




1 

1 
















i 

1 




1 if 




1 

1 






| 










1 




1 if 




1 






1 










1 




1 if 




1 






1 










1 




1 if 




1 






1 










1 




1 if 




1 






1 










1 




1 if 




1 






1 










1 




1 if 




1 






1 










1 




1 ■* 




1 






1 










1 




1 if 




1 






1 










1 




1 if 




1 






1 










1 


X 


1 if 




1 
















1 


•— » 


1 -if 




1 






LU 










1 


(X 


1 if 




1 






UO 










1 


1 — 


1 -if 




1 






o 








O 


1 


<r 


1 if 




1 






a 








LU 


1 


X 


1 if 




1 






CC 








a 


1 




1 if 




1 


XX 


XX 


o 








UJ 


1 


o 


1 -if 




1 


CL 2 


CL "2 


CL 








Ui 


1 


2 


1 if 




1 




~r ZZ2 








— * 


— * 


1 


*— » 


1 if 




1 


UJ - 


UJ - 


o 










1 


1 — 


1 if 




1 


K< 


H3 


TF 


— * 






a. 


1 


X 


1 if 




1 


-zs 




>— 4 


— o 




o 




1 


o 


1 if 




1 


<TG 


30 


X 


X—* 




— 




1 


•—* 


1 if 




1 




3:uc 


o 


OC * 




— 


LU 


1 


LU 


1 if 






OX 


OX 


LU 


OX 




m 


OC 


1 


2 


1 if 




X 


or - 


(X - 


X 


XcL 




— 


OJ 


1 




1 if 




H-4 


x< 


xo 


o 


-O 




o 


X 


1 


LU 


1 if 




UL 


-o 


-o 




(3,2 




•> 




1 


f- 


1 if 




h- 


uro 


30 


X 


— X - 




— • 


CO 


1 


< 


1 if 




<x 


oo 


30 


o 


30 1 — 




cm 


LU 


1 


h~ 


1 if 


o 


2T 


OX 


OX 


LU 


30 LU 




«■— 


a 


1 


OO 


1 if 


o 




o - 


o - 




OCLO 




o 


o 


1 




1 if 


3 


J 


x< 


xo 


X 


o — - 




► 




1 


LI. 1 


1 if 


o 




- 2 




"2 


X^l * 




— * 


UJ o 


1 


X 


1 if 


X — 


X 


30 


30 


OcO 


->>^ 




r-«i 


X ca 


1 


I — 


1 if 


-O 


Z) 


enc 


lu: 


Ui 


c?cl a. c\i 




— 


»— i | ■ 


1 




1 if 


oo 




ox 


ox 


LLX 


2 1 — * 1 — ' 




OO 


h- X 


1 


1— 


1 * 


id 


X 


cr - 


a. - 


o^h 


O — — H - 




•*.30 


o o 


1 


o 


1 if 


Oo 


o 


XX 


XX 


h— 


a: 33m 


2 | 


- -r 


O X 


1 


a. 


1 >f 


OCX 


•-H 


-CL 


-CL 


f-O 


X 3 30 


3 1 


f— 


a: - 


1 


i — 


1 if 


2 - 


r- 


— X 




XO 


-O O — ' 


1 — 1 


XO 




1 


o 


1 if 


-3 


< 


— MLJ 


NU.) 


<Xc£ 


X003 


UJ | 


1 — ‘I — 


or 


1 


o 


1 >f 


Ol X 


?: 


r>M— 


Oh- 


y 


CLXX< 


a i 


<T 


o cn 


1 




1 if 


CLO 


cc 








X ~ -UJ 


•* I 


ecu 


►-H 3 


1 


m* 


1 if 


XCl 


o 


mx»— * 


X *— • 


XO 


U J O O CL 


« | 


1— ' 0 


h- Cl 


1 




1 if 


LU 2 


X 


X^r_J 


"5N 


ql <r 


1 ii 


2 | 


(X) 


< h- 


1 




1 if 


h- - 




f- ON 


30 


LUCL 


- - -ir 


3 1 


~T ,— . 


k* 2 : 


1 


o 


1 if 


-X 


< 


^ w 




♦ — X 


xxx • 


f— 1 


u o 


o a 


1 


►—4 


1 if 


x a 


QC 






LUj— » 


2 cl a . o 


- UJ 1 


o • 


CL X 


1 


h- 


1 if 


zz tn 


1 — 


HJH 


Jh 


Q 3 


oxx-h 


O 1 


• f- 




1 


< 


1 if 


<— 3ui 




• XX 


XX 




— UJUJ-Ji 


— — 1 


o 


U 


1 


N 


1 >f 


u-h- 


UJ 


O) — f— 


h- X 


OO 


i — p— 




— * • 


o — 


1 


>— • 


1 if 


cc<r-^ 


O0 


UJ<< 


<t<r 


XX 


O- 


| 


if 3 


if 


1 


21 


1 if 


com li j 


ex 




XX 


*— • 1— • 


UQ-- 


if if 1 


-OCO 


on 


1 


►-H 


1 if 


-nr K- 


UJ 


XOU) 


oo 


U.U) 


UJLjOh- 


— — 1 


or-' sjui 


4£L(NJU- s U 


1 


h- 


1 if 


>0-1'^ 








3 


3 LU LU u J *U vO 1 


X •"* —> 


HI <—4 ) — » 


1 


CL 


1 if 


^CLIX 




UJJ 


30 




0.000 




o x 


—• z: 


1 


a 


1 if 


O < 




UJJ 


-J3LL 1 


OOO || 




ujxO— < 


QQmLU 


1 




1 if 


LU 


1 


— •<< 


<< 


►— • | 


2 UJ LJU | 


1 ’1-H 


CM-l” 


1 


X 


1 if 


h-IJ 


1 


-UOOJUU 


1 


3330 




»-•« X 


or x~ 


1 


C) 


1 if 


-.JJ 


1 


oo 


O 1 


3331- 


« H-H | 


r/LLOO 


'DJOrr 


1 




1 if 


ex <r <f 


1 


X 


3 


X 1 


<r O T ! j or r/ i 


X — »0 O 


oo 2 


1 


Q 


1 if 


.20 0 


1 


*—• 1 


U» 1 


Ui I 


0003-2 2 1 




UJ 


1 


X 


1 if 




1 






1 








X 


1 


LU 


1 if 




1 






1 








r— 


1 




1 if 




1 






1 






O O 


o 


1 




I ;f 




1 






1 






0 n 


CO 


1 




i if 

















OO'J’ rooo iU'JUUOOuOO UUO OOO UUUJ 



131 



CALL CPLRFA(3WX»Qir. , NROWQ* NCGLQ ) 

U. ITE U,*) • END OF POLE PLACEMENT ROUTINE • 

aRITE ' ' 

WRITE ( 6 * * ) • THE C REQUIRED FOR THIS ASSIGNMENT IS 
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CALL MMUL ( f ,C , .vIRQU F , NCGL F , NC 0 LC , F C i CGN02L60 

CALL M . v L L (B.FCtNRCWB.NCJLE.NCOLC.BFCJ CC.\'03c70 

DO 180 1 = 1 t NRO W A CCNC3L60 

DO 170 J=1,NC0LA CCNC3690 
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